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Rational functions with maximal radius of absolute monotonicity 

Lajos Loczi and David I. Ketcheson 
Abstract 

We study the radius of absolute monotonicity R of rational functions with numerator and denominator of degree s 
that approximate the exponential function to order p. Such functions arise in the application of implicit s-stage, order 
p Runge-Kutta methods for initial value problems, and the radius of absolute monotonicity governs the numerical 
preservation of properties like positivity and maximum-norm contractivity. We construct a function with p — 2 and 
R > 2s, disproving a conjecture of van de Griend and Kraaijevanger. We determine the maximum attainable radius for 
functions in several one-parameter families of rational functions. Moreover, we prove earlier conjectured optimal radii 



o 



in some families with 2 or 3 parameters via uniqueness arguments for systems of polynomial inequalities. Our results 
also prove the optimality of some strong stability preserving implicit and singly diagonally implicit Runge-Kutta 
methods. Whereas previous results in this area were primarily numerical, we give all constants as exact algebraic 
numbers. 

(N 

1. Introduction and aims 

A smooth function -0 : K — J' K is said to be absolutely monotonic at a point if ijj and all of its derivatives 

r*^ are non-negative there. The radius of absolute monotonicity R{ip) G [0, +oo\ is defined as 

> 

^ R{tj-') — sup ({r G [0, +oo) : tJj is absolutely monotonic at each point of [— r, 0]} U {0}) . 

B 

The radius of absolute monotonicity of polynomials and rational functions plays an important role in the 
analysis of positivity, monotonicity, and contractivity properties of numerical methods for initial value 
" ' problems and is often referred to as the threshold factor in this context [l] [oj |12[ Specifically, 

^ I the maximal positive or contractive step-size is given by i?(?/')/io, where "0 is the stability function of the 

numerical method and Hq is the maximum step-size under which the corresponding property holds for the 
\^ explicit Euler method. 

It is therefore natural to consider the problem of finding a function ijj that achieves the maximal radius of 
absolute monotonicity within a given class. In this work, we study absolute monotonicity of rational functions 
that correspond to the stability functions of certain implicit or singly diagonally implicit Runge-Kutta 
methods. 

\ I A Runge-Kutta (RK) method of s e N+ stages is defined by its coefficients, an s x s matrix A and an 

I> s X 1 vector b [2] . The stability function of the method is 

X ^ ,Ab/ N det{I - zA + zlb^) , ^, ,^ ^, 

where 1 is a column vector of length s with all unit entries and / is the s x s identity matrix [5^ . The order 
of a method, denoted by p G N+, indicates how accurately the computed solution approximates the exact 
solution in an asymptotic sense. The stability function of an RK method of order p must approximate the 
exponential function to at least order p near the origin. 

In an implicit Runge-Kutta method (IRK), all entries of A may be non-zero. An important subclass of 
the IRK methods are the singly diagonally implicit RK (SDIRK) methods, with A lower triangular with 
identical diagonal entries. For explicit RK methods, A is strictly lower triangular. 
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We now define the classes of rational functions to be studied. For m, n S N, let Hm denote the set of real 
polynomials of degree at most m: 

n™ = I £ ft-^' : f3, eR, j = 0,1,... ,m\, 

and let Tl„i/n,p denote the set of all real (m, n)-rational functions approximating the exponential to order p 
near the origin: 

P 

Let Tlm/n,p denote the elements of Tlrn/n,p whose denominator has (at most) a single, non-zero real root: 



n,„/„,p = {V : = ^, P e n„, ^ Q e n„, ^{z) - exp(z) = 0(2^+^) as z ^ 0}. 



^ r ^y^j 1 

n,„/„,p = jV' e n,^/„_p : i/'(z) = _ p e n„,a e m|. 

Obviously, for every Q < m < rh, Q < n <n and 1 < P < P , we have n„j/„ p C l^m/n,pi l^m/n,p C Tlfh/n,p and 
If A and 6 correspond to an IRK method of s stages and order p, then 

e n,/,,p, (1.2) 

while for SDIRK methods, 

V-^'" e n,/,,p. (1.3) 

For explicit RK methods we have ■0'^'^ e Hg. 

A thorough study of polynomial approximations to the exponential with maximal radius of absolute 
monotonicity can be found in |9j [g]. We are interested in determining the maximal radius of absolute 
monotonicity that can be achieved among the stability functions of IRK or SDIRK methods of a given order 
p. Therefore, for non-empty sets Ilm/n,p and Tlm/n,p, let us define the quantities 

Rm/n^p SUp{i?(-0) : tp e U^/n,p} 
Rm/n,p SUp{i?(-0) : i) e U^/n,p}- 

We will focus on the cases for which m = n = s. 

The seminal work on this topic is [12| , in which an algorithm is presented for computing the radius of 
absolute monotonicity of a rational function, and many properties of the radius of absolute monotonicity 
are proved. The determination of R{ip) is not trivial even for a single rational function Tp, so the difficulty 
of obtaining Rm/n,p or Rm/n,p ioi a particular (rn,n,p) triple ranges from fairly challenging to currently 
impossible. Nevertheless, some patterns in numerically computed values have led to important conjectures. 

In order to be able to fully describe these conjectures, we recall the concept of radius of absolute 



monotonicity of a Runge-Kutta method 10 , denoted by R{A,b). The quantity R{A,b) is also referred 
to as Kraaijevanger's coefficient [s], or SSP coefRcient For a RK method with coefficients A,b, define 

if e]R(- + l)x(s-M) by 

K = K{A,b):= ^ [] 

Let 1 denote now the vector (1,1,...,1)^ e M'*+^. Then the radius of absolute monotonicity of the RK 
method is 

R{A,b) := sup{r e M : Vp G [0,r] 3{I + gK)'^ , gK{I + gKy^ > and gK{I + gK)'^! < 1}, 

where vector and matrix inequalities are understood componentwise. Notice that R{A, b) > 0. Absolute 
monotonicity of a Runge-Kutta method implies absolute monotonicity of its stability function |10j ; thus we 
have 

R{A,b) < RitP^^"). (1.4) 

The coefRcient R{A, b) plays the same role in numerical preservation of positivity and contractivity for 
non-linear problems that the coefficient R{^^'^) plays for linear problems [4j. 
The following conjectures served as motivation for our work. 
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Conjecture 1.1 ([l2| p. 421]). For m,n eN~^, Rrn/n.2 = fn 



In 12 , this conjecture was proved for all to > 1 with n = 1 or 2. In the special case m ^ n = s E N+, the 



conjecture claims 



/s,2 



2s. 



(1.5) 



For A,b corresponding to an RK method with s stages and order p, we have R{ip ' ) < Rs/s,p- Moreover, 



R 



s/ s,p 



is a non-increasing function of p for fixed s; therefore (1.4) and (1.5) together would imply that 



R{A, b) < 2s for all Runge-Kutta methods that are more than first order accurate. Indeed, evidence in the 
literature reinforces belief in the bound ( 1.5 ); we have the following conjectures based on numerical searches. 



Conjecture 1.2 ([sj Conjecture 3.1]). Let A,b denote the coefEcients of an SDIRK method of order 
p>2 with s > 1 stages. Then R{A, b) < 2s. 

Conjecture 1.3 ([s]). Let A,b denote the coefEcients of an imphcit RK method of order p>2 with 
s > 1 stages. Then R{A,b) < 2s. 



It is thus very surprising that, as we will see. Conjecture 1.1 — and thus equality (1.5) — does not hold 
for m = n = 3. An immediate consequence is that Conjecture 1.3 cannot be proved by analyzing Rs/s.2- A 
weaker conjecture that arises naturally is 

Conjecture 1.4. For each s > 1, Rs/s,2 — 2s. 



Conjectures 1 1.2 1.3 and 1.4 were previously proved only in the cases s = 1 or s = 2. 

We have presented these conjectures in the order they were formulated. By also taking into account that 

• Rs/s,p < Rs/s,2 < Rs/s,2 [p > 2), 



R{'ip^'°) < Rs/s,p for A, b corresponding to an SDIRK method, 



• there exists an SDIRK method (consisting of the repetitions of the implicit midpoint method) with 
R{A,b) — 2s (^ formula (3.1)]) whose stability function 

{'^ 2s) 

satisfies R{^p^'^) — 2s, we have the following implications for each value of s: 

Conjecture |l.l| with m = n = s Conjecture |1.4| Conjecture |1.2| 

and 

Conjecture 1 1.1 1 with m = n = s Conjecture |1.3| Conjecture |1.2 [ 

As for the p — 3 case, we pose the following new conjecture, which is stronger than |3, Conjecture 3.2]. 

Conjecture 1.5. For each s >2, Rs/s,3 = s — 1 + Vs^ — 1. 

The present work is devoted to determining the exact values of Rs/s,p and Rg/s.p for certain 2 < s < 4 
and 2 <p <1 pairs. Our results lend some support to Conjectures |1.2| and |1.4[ since we prove each of them 



for the cases s = 3 and s = 4. We also prove Conjecture 1.5 for s < 4 (simultaneously proving the 3 < s < 4 
special cases of [3| Conjecture 3.2]). 

The structure of the paper is as follows. 

In Section [l.l I we recall some theorems from 12 10 on which our computations are based. 



All of our new results are overviewed and interpreted in Section [2j the proofs are given only in Sections 
3]-[9] Section [2] has been written so that the casual reader need not refer to any later sections. In Section 
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2.5 we give some auxiliary, but, in our opinion, related and interesting results that we could not (yet) tie 



to the main pieces of the puzzle (e.g., to Conjecture 1.4), along with a few remarks about the successful (or 



failed) proof attempts and techniques. In Section 2.6 we introduce some compact notation for the algebraic 
numbers that will be prevalent in this work. 

The interested reader will find the actual proofs of our theorems in Sections |3]-[9j Generally, we proceed 
from smaller s values to larger ones, and if s is fixed, from larger p values to smaller ones, i.e., sections 
are ordered roughly in increasing difficulty within both the Hs/s,p and in the Tls/s,p classes. These sections 
represent the fruits of several dozens of pages of computations — or of a few hundred pages, depending on the 
level of detail — so we had inevitably to omit some parts of some proofs. In Section |9j a few more algebraic 
expressions, mentioned only implicitly in the proofs, are collected to enable the reproducibility of certain 
longer computations. 

Let us close this introduction with a remark explaining why the p = 1 case is exceptional. 



Remark 1.6 (on the p = 1 case). For ip G Tlm/n,p with p >2, we have R{ip) < +oo 11 '. On the other 
hand, for ip{z) — 1/(1 — z), we have ip S IIo/i.i and R{ip) — oo. One may ask whether there exists ip G 
lig/s,! such that i?('0) — oo and ip ^ Il(^s-iy^s-i),i (a-nd correspondingly for the sets Ils/s,i)- The answer is 
atRrmative. For any s G define 

1 1 



2s 2s(l - 2zy 



and for any 2 < s G N, consider 



1 4^ - 1 ^ 21-'" - z' 

m— 1 



It is easily seen that both ipg and ips can be represented as rational functions with numerator and 
denominator having degree exactly s. Moreover, -ipsiO) — i^'si^) = 1 = ^s(O) ~ V's(O); so -tps G tts/s,i 
tps G ^s/s,i \ ^s/s,i (of course, ipi € n^/i ^ = IIi/i ^j. Direct differentiation shows that for all k G N and 
a; < we have ■ipi''\x) > and ^i''\x) > 0, hence R{'ijjs) = R{4's) = +oo. Consequently, we are going to 
consider only the p > 2 case in this work. 



1.1. Some general results on the radius of absolute monotonicity 



Let us briefly summarize some useful theorems from 12 10 that will frequently be used in this work. In 



12 , the following assumptions are made on the rational function i/j: 



® V = where P G Um and Q € n„ with some m,n gN, but ijj is not a polynomial, 
© P and Q have no common roots, 
® P(0) = Q(0) = 1. 

Remark 1.7. We will see that removing these assumptions (i.e., not excluding removable singularities a 
priori, or considering, when setting up the form of the families of rational functions, the case P(0) = Q{0) = 
as well with interpreting V'(O) = 1 as limp ^ = 1 in the order conditions) does not make a difference in the 
optimal values Rs/s,p and Rs/s,p in the classes we are going to consider. Nevertheless, these assumptions are 
convenient in allowing us to immediately use results from [12 J. 



Definition 1.8 ( 12 Definition 3.2]). Suppose that ip satisfies Assumptions ®— (D above. Let A+(V') 



denote the set of poles of ip with non-negative imaginary part. If a G A^{ip), we set 
I{a) :— {x gM. : a G A'^ is the unique pole closest to x}. 
The disjoint union of these intervals is the set R with only finitely many exceptions. Now we let 



— inf /(ao) if G /(ag) for some positive real pole ao G A'^(^p) n (0, +oo), 
otherwise. 
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Note that B{ip) is determined solely by the location of the poles of ip. The significance of this quantity is 
highlighted by the following theorem. 



Theorem 1.9 ( 12 Theorem 3.3]). Suppose that ip satisfies Assumptions ®— (D above. Then 



Theorem 1.10 ( [12[ Corollary 3.4]). Suppose that tp satisRes Assumptions ®— (D and ip has no positive 
real pole. Then R{ip) = 0. 

Additionally, the negative real roots of the derivatives of "0 form upper bounds on Riip), as shown by the 
next theorem. 

Theorem 1.11 (jl2j Lemma 4.5]). Suppose that ip satisfies Assumptions ®—®, and tp'^^^ {x) = for some 
> a; e K and £ e N. TJien Riip) < -x. 

We will also use the fact that under certain assumptions absolute monotonicity at the left endpoint of an 
interval implies absolute monotonicity on the whole interval. 



Theorem 1.12 ( 
are polynomials anc 



10 



Lemma 3.1]). Let ip — ^ be absolutely monotonia at some x < 0, where P and Q 



"Q^has no zeros in (x, 0]. Then R{tp) > —x. 



Remark 1.13. In ^12^, formula (4.3) introduces an auxiliary quantity 

L{x) — max (0, m — n + 1, max{. . .}) 

that is used in their algorithm to compute R{'p) for a given rational function. It may happen however that 
the {. . .} set above is empty (corresponding to the SDIRK case, for example), when a correct interpretation 
of this max0 is — oo (or, say, 0). 



2. Main results 

In this section we state the main results and describe the approach of the proofs. The proofs themselves 
are deferred to Sections [3l-|9l 

For a given set Ils/s,p or Tls/s,pj let q = q{s,p) denote the number of parameters needed to describe the 
set. Generically, the class Tls/s,p can be written as a family in q = 2s — p parameters, while the class Tlg/s.p 
can be written as a family ing = s + l— p parameters. In the cases where q = 0, these sets contain a finite 
number of members ("finitely many" is of course understood in the sense of functions, i.e., a normalized 
representation is chosen: if ip — ^ with V'(O) — 1, then P(0) = Q(0) = 1). Our investigation is restricted to 
sets consisting of g-parameter families of rational functions with q = (3 such cases), g = 1 (5 cases), q = 2 
(2 cases), q = 3 {1 case), or g = 4 (a lower bound in 1 case). 



2.1. A lower bound on R^/ 



3,2 



1.1 



In Section 4.3 we construct a function ip e 113/3 2 with i?('0) > 6.7783 > 6. This shows that Conjecture 
does not hold for m — n = S. However, the exact value of i?3/3,2 is still unknown, because in order to describe 
all rational functions ip in 113/3 2, we need 4 parameters, rendering the (exact or numerical) optimization 



within this class impossible with our current techniques. In contrast. Theorem |2.7| in Section [2. 4| asserts that 
-^3/3,2 = 6. 

2.2. Determination of Rs/s,p for s < A, q < 1 
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2.2.1. Known values of and R2/2,p For p > 3, Hi/i.p is empty. The set II1/1.2 consists only of the 

function z jz^r/^ with i? = 2, so -Ri/1,2 — 2- 

For p > 5, n2/2,p is empty. The set 112/2,4 consists of a single function: the (2, 2) Pade approximation of 
the exponential function. This function has '0'®-'(O) = 0, so by Theorem 1.11 we have i?2/2,4 = (but see 



Remarks 2.1 and 2.3 as well). The same result was deduced in a different way in [12, Theorem 5.1]. 

If ip = e 112/2, 2j E^nd P and Q have no common roots, then jl2^ Section 6.3] says that R{ip) < 4. If P 
and Q have a common root, then it is easily seen that ip g ^1/1,2 also, hence R{ip) < 2. For the function 

^(^) = {l^W have = 4, so i?2/2,2 = 4. 

The exact value of R2/2.3 'was previously unknown; it is determined in Section |3j thus completing all cases 
with 1 < s < 2. 



2.2.2. JVew exact values of Rs/s,p For the cases in which p = 2s, the set Ils/s,p has only one member: 
the (s, s) Pade approximation to the exponential. The value of Rs/s,p can be obtained by computing R{ip) 
for this unique member; see also Remark |2.3[ 



For the cases in which p — 2s — 1, the set ^s/s.p is a one-parameter family. In these cases, the method of 
proof we use is the following procedure (possibly in an iterative manner). 

1. Conjecture an optimal parameter value a* by inspection of B{ipa) and the first several derivatives of 

■0a- 

2. Rigorously exclude all the parameter values E\{a*} by appealing to Theorem |1.11| and the 
intermediate value theorem. 

>0, 



3. 



Explicitly compute a formula for the fc'^ derivative (/c € N arbitrary) and prove that 1/;^? 

where x* > is as large as possible due to Step 1. 
Decimal approximations of proven optimal values are given in Table [2?T| along with some other properties. 
Exact values are given in the sections indicated. 



Remark 2.1. We remark that the optimal ip functions corresponding to the {s,p) = (2,2) and (s,p) 



(2,3) rows of Table 
and so R < B. As 



2.1 



arc also elements of the corresponding class Hg/s.p- Hence for those functions B — +00, 
for the optimal ■(/' corresponding to {s,p) — (2, 4), we have checked that for < k < 100, 
^('=)(0) = if and only if k is divisible by 6. 



Remark 2.2. 
TaWe 6.1]. 



Our computation of i?2/2,3 confirms the corresponding numerical result presented in [12 



Remark 2.3. The classes Tls/s,p where p = 2s each contain a single element: the {s,s) Pade approxi- 



mation to the exponential. For s even, it was shown that Rs/s.2s — in 12 Theorem 5.1]. For s — 3, the 



(s,p) 



R 



s/ s.p 



deg(a*) <ieg{R^/^^p) 



Section(s) 



Obstacle (with R := Rs/s.p) 



(2,2) 
(2,3) 
(2,4) 
(3,5) 
(3,6) 
(4,7) 
(4,8) 



:4t 

: 2.732050* 
:0t 

: 2.301322 
: 2.207606* 
: 2.743911 

:0t 



30 



1 
2 
1 
6 
3 

30 
1 



2.2.1 



T7I 



IT 

5 land 



emark 



il,{-R) = ip'i-R) = 
tP'{-R) = 

_R = B = or i/)(6'")(-i?) 

R = B 

R = B 

R = B 

R = B = 



Table 2.1. Optimal Rs/s,p values, together with the algebraic degree of the optimal parameter choice 
within the given parametrization, the algebraic degree of i?s/s,p> the section number in which Rs/s,p is 
given as an exact algebraic number, and the factor that limits the optimal value (see Section 1.1). 



Superscripts ^ and ^ indicate, respectively, that the optimal Rs/s,p value was already proved earlier exactly 
or was presented earlier numerically. A — means that the corresponding Tls/s,p class contains only finitely 
many elements hence the parameter a is not present, whereas a + symbol denotes that the class cannot he 
described by only one parameter. As for the integer m in the {s,p) — (2, 4) case, see Remark 2.1 
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optimal value i?3/3.6 ~ 2.2076 is presented numerically in 12 Table 5.3], and in Section 4.1 we have given 
its exact value as an algebraic number of degree 3. 

2.3. Determination of Rs/s.p for s < 4 



The known and newly obtained Rg/s.p values are summarized in Table 2.2 Now let us briefly interpret 
our results in these classes. 

For 1 < s < 4, the classes Hg/s.p are em pty when p > s + 2; the trivial class is mentioned in Section 



2.4 



So by taking into account Table 2.2 we have a complete description of the optimal Rg/s.p values for 



s < 4 and P > 2. 

The classes 113/3 2 £^nd II4/4.2 will be treated separately in Section 2.4 so as to be able to give a detailed 
exposition of their (more challenging) proofs. 

The Rs/s,3 values for 2 < s < 4 support Conjecture [l.5[ being stronger than [3| Conjecture 3.2]. In other 
words, we have proved |3j Conjecture 3.2] for 3 < s < 



Remark 2.4. On one hand, according to (1.4), if A and b correspond to any SDIRK method of s stages 
and order p, then we have R{A,b) < Rs/s,p- On the other hand, an SDIRK method is constructed in |3| 
for each pair {s,p) satisfying 2 < s < 4 and 2 < p < 3 such that R{A, b) — Rs/s,p- Based on this, one might 
suspect that Rs/s,p J's equal to the optimal R{A, b) radius for s-stage, order p SDIRK methods for each 
s,p > 2. However, Remark \2. 5\ shows that this is not the case in general. 

Remark 2.5. There exist exactly 3 functions in the set 113/3 4. These 3 rational functions are also 
mentioned in |51 Section 3.4.1]. The method satisfying the non-negativity condition K > in Jsl is the one 
whose stability function yields the optimal i?3/3,4 value. This function has radius of absolute monotonicity 
-^3/3,4 ~ 3.2872, but the corresponding optimal SDIRK method has only R{A, b) 1.7587. 

Remark 2.6. The optimal ipa* in the 1X4/4 4 class is different from the stability function of the optimal 
method of ^ Section 3.4.2] in this class obtained by numerical search. In other words, according to the 
numerical tests in the optimal R{A, b) in the SDIRK s — p — A case is w 4.2081 < i?4/4,4 (c.f Remark 



2.5) 



Finally, we give some additional information on the structure of the optimal ip functions in Table |2.3| 

2.4. Determination of Rs/s.2 for 3 < s < 4 

The case p = 2 is of particular interest in light of the conjectures presented in Section [T] 



(«.p) 


D 


deg(a*) 


deg(i?a/s,p) 


Section(s) 


Derivatives that vanish at —Rs/s,p 


(2,2) 


= 4t 


1 


1 


I2.2.I 




i = 0,l 


(3,2) 


= 6 


+ 


1 








e = 0,1,2 


(4,2) 


= 8 


+ 


1 


7 anaT9.2 


^ = 0, 1,2,3 


(2,3) 


2.732050* 




2 


Remark 


2.1 


e = i 


(3,3) 


^ 4.828427 


2 


2 




6.2 




i=l,2 


(4,3) 


^ 6.872983 


+ 


2 




8.3 




^ = 1,2,3 


(3,4) 


^ 3.287278 




9 




6.1 




e = o 


(4,4) 


^ 5.167265 


9 


9 




8.2 




1 = 0,1 


(4,5) 


^ 3.743299 




12 




8.1 




1=1 


Table 2.2. 


Optimal Rs/s,p 


values. 


together with the algebraic 




iree of the optimal parameter value in 



the given parametrization and the algebraic degree of Rg/s.pj section number in which they are gi ven as 
exact algebraic numbers, and the first few derivatives of ^ that vanish at —Rs /s.v (see Theorem 
symbols ,+,— have the same meaning as in Table 2.1 



1.11). The 
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(2,2) 
(3,2) 
(4,2) 

(2,3) 

(3,3) 

(4,3) 
(4,4) 



Form of the optimal ip or its derivative 



V(^) = (1 + ^/4)2/(1 -2/4)2 
V(2) = (l + 2/6)-V(l- 2/6)3 
V>(2) = (1 + 2/8)V(l - 2/8)4 

i,'{z) = (l + I [V2^ - l) ^) / (l + i [V2^ -3)zy 
\ + \ {VW^ - 2) ^) ' / (l + i (v^3^ - 4) z) 
1 /V4^-3) .)'/(l+i^(V4 



V'(2) 
l/''(2) 



1-52 



,/,(2) = (1 + 0.1935262)2 (1 + 0.2244162 + 0.043763822) /(I - 0.09713312)^ 



Table 2.3. Each row shows the appropriate derivative of the optimal element in Hg/s.p having multiple 
roots. In the {s,p) = (4,4) case, approximate constants are used. In the {s,p) S {(3,4), (4,5)} cases, no 
multiple roots were found among the Grst few derivatives of the optimal ip functions. 



Theorem 2.7. Fix 1 < s < 4. Then we have Rs/s,2 — 2s, and the unique tp G Tls/s.2 that attains -R(^) 
2 s is 



(2.1) 



The values i?i/i.2 = 2 and R 



?2/2,2 = 4 were determined already in |12| and formed part of the basis for 

that the optimal elements in Hj^/i 2 112/2.2 si's 

R2/22 = 4. For the sake 



2.2.1 



Conjecture |1.1| Moreover, we have seen in Section 

also elements of IIi/i 2 and 112/2,2, respectively, so i?i/i,2 = ^1/1,2 = 2 and i?2/2,2 = ^2/2,2 
of completeness, we give two short and direct proofs of the equality R2/2.2 = 4 at the beginning of Section 
|6j illustrating the two proof strategies we use in this work, and simultaneously proving the s = 2 case of 
Theorem 1 2. 7| The 3 < s < 4 cases are new. The proof for s = 3 is given in Section [STS} while the s = 4 case 
is described in Sections [7| and [9^ Theorem |2. 7| confirms Conjecture |1 .41 for s < 4, furthermore, it also proves 
Conjecture 3.1 in [sj Section 3.2] for 3 < s < 4 (the truth of this conjecture for 1 < s < 2 being already 
established in [s]) regarding SDIRK methods that are optimal with respect to R{A,b). 

We give here just the common ingredients used in both proofs (s = 3 and s = 4), then the proofs are 
finished in the corresponding sections as indicated above. 



For s > 3, any ip € ^s/s.2 can be written as 

1 + (1 - a{l)) z+{\- a{l) + a^Q) + 



(1 - az)" 



(2.2) 



with some a, 03, 04, . . . , e M. Let us denote the numerator of ( 2.2 1 by P{z) and set Q{z) 1 — az, so that 



The next step is to exclude the non-positive a values. Since Theorem 2.7 will be proved via the uniqueness 
argument 



Riip) = 2s implies the unique form (2.1 ) of tp, 
and 2s > s — 1, the assumption of the lemma below is justified. 



Lemma 2.8. Let Tp S Ils/s ,2 be given in the form (2.2), and suppose that R{ip) > s — 1. Then the 
parameter a appearing in (2.2) satisfies a > 0. 



Proof. If a — 0, then ^ is a polynomial approximating the exponential function to order p — 2 near the 
origin. For such polynomials, R{ip) < s — 1 [9] Theorem 2.1]. For a < 0, tp has no positive real poles, so we 
will apply Theorem 1.10 to show that R{ip) = 0. We need only to verify Assumptions ®-@. Assumptions 
® and ® are automatically satisfied by functions of the form (2.2), so we are done if Assumption ® is also 
fulfilled. If not, then except at z = 1/a, tp can be expressed as P/Q, where P and Q have degree 's < s and 
no common roots. Clearly, s > 0, since p — 2 still hold s. Mo reover, P(0)/Q(0) ~ P{Q)/Q{0) — 1, so we can 
assume P{0) — Q{0) = 1. Now application of Theorem 1.10 to P/Q shows that a < implies R{tp) =0. □ 
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Now a useful necessary condition is derived based on the fact that if a polynomial in k is non-negative for 
all k G N, then its leading coefficient is also non-negative. Assumption a > in the next lemma is guaranteed 
of course by the previous lemma. 



Lemma 2.9. Suppose that ip = given by (2.2) with a > is absolutely monotonic at some a:: < 0. 
Then 



1 ■^-^ 1 /s ^ I + k — m\ f 1 — ax 
hm ^-^ 2^ 



fc->-oo k^ 



m=0 



m 



s - 1 



a 



(2.3) 



Proof. From (2.2) it can be shown by induction that the k derivative of t/; (fc G N) is 



V I I Qs+k ^! 

^ m=0 



s-l 



Q' 



p{m) ^ 



(2.4) 



Suppose that tp = is absolutely monotonic at some a; < 0. We have now —Q'{x) = a > and Q{x) > 0, 
so (2.4) implies 

min(fc,s) 



lim — 

fe— >-oo k 



- y - 



m— 



1 + fc — m 
s~ 1 



Q'ix) 



P(")(a;) > 0. 



(2.5) 



Note that the above limit always exists since the sum is a polynomial in k of degree at most s — 1 and the 
sum is finite. □ 



After these preparations, we show in Section 6.3 for s = 3, and in Sections [7] and 9.2 for s — 4, that 



(2.1) is the unique solution of the following non- linear system of polynomials in the variables a > and 
as e M: 

^(k) (_2s) > for = 0, 1, . . . s - 1 (2.6) 



and 



(2.7) 

The strength of the above uniqueness-type argument is that we could handle more than one parameter 
within this framework. The disadvantage is, however, that the optimal value (here R = 2s) must be known 
in advance. The reason for separating indices < A: < s — 1 and fc > s in the above non-linear system will 



be explained in Remark 2.23 



2.5. Further results and questions 

In this section we present some intermediate results that may prove useful in future studies of the radius 
of absolute monotonicity. 



The motivation for Section 2.5.1 came from Conjecture|1.3| Let A, b denote the coefficients of an arbitrary 



RK method; then i?(A, 6) > =^ A>0 |4, Observation 5.2] (again, matrix inequalities are understood 



componentwise). The trace inequality (2.10) for non- negative matrices is used to derive a certain non-linear 
relation (2.11 ) between the first few coefficients of the polynomial appearing in the numerator of the stability 
function ( |1.1[ ). 

Next we turn our attention to Conjecture |1.4[ Let us fix s > 3 and p = 2, choose an arbitrary S DIRK 
method with coefficients A,b, and consider its stability function -0 g ^s/s,2 as described by formula (2.2 1. 
Let ao,ai and a2 denote the first few coefficients of the numerat or of ip, that is oq := 1, ai := 1 — a(^) and 
a2 :— \ — +0^(2), then a simple computation shows that (2.11) is satisfied — in fact, with equality — 
despite the fact that the non-negativity condition A>Q here is not assumed. If 3 < s < 8, then Lemma[2T6] 
in Section [23^ gives a remarkable uniqueness result for polynomials: the unique polynomial given by Lemma 



2.16| is identical to the numerator of the conjectured optimal and unique function appearing in Theorem 
2.7| (or in Conjecture 1.4). Apart from the fact that we could prove Lemma 2.16 only for 3 < s < 8, the 



missing link is the following: the uniqueness result in Theorem 2.7 is essentially obtained under the condition 

^(fe)(_25)>0 for fceN, (2.8) 



Page 10 of [42 



LAJOS LOCZI AND DAVID I. KETCHESON 



whereas in Lemma 12.161 we assumed 



p(fc)(_2s)>0 for fc = 0,l,...,s. 



(2.9) 



We are very curious whether a result similar to Lemma [2.161 could lead to a proof of Theorem ] 2. 7| for general 
s values, tha t is. to a proof of Conjecture |1.4[ 



In Section 
Section 12.5.4 



2.5.3 



we give a formula for the elements of the Iis/s,p class for general s and p values, while in 



wc briefly touch on some of our proof techniques and their possible extensions. 



Remark 2.10. We have found the foUowing conjecture whose assum ptio n — with pm = P^"^\x)/m\, q = 
Q{x) and q — Q'{x) = —a for some x E [—2s, 0] — is similar to the sum in {2.4) apart from the factor (— 1)*^^™ 
and the fact that the Pm quantities can be independent of one another. The conjecture has been proved by 
Mathematica for 1 < s < 10, and we have a "manual" proof for s — 2. We do not know whether the conjecture 
could be used in a transition from conditions \2.8) to (2.9). 



Conjecture 2.11. Fix any s G'N~^,0^q€R, 0<q£M. andpm € M ('m = 0, 1, . . . , s + Ij, and suppose 
that for each fc = 0, 1, . . . , s + 1 we have 

min(fc,s) , 1,1, \ 



m=0 



Tien pq>Q and for to = 1, 2, . . . , s + 1, we have Pm = Po (^) ' ( | ) • Moreover, tiie above sum is equal to 
Po, if k = 0, and 0, if k = 1, 2, . . . , s + 1. 



2.5.1. Trace inequalities for non-negative matrices and the numerator of the stability function Conjecture 



1.3 in the s — 3, p — 2 case claims that R{A, b) < 6 for any RK method with coefficients A, b. The following 



general lemma was discovered while investigating this conjectured bound. Due to the remark in the beginning 
of Section [2?5| A> will be assumed throughout the current Section [2.5. 1[ 

Let us fix a positive integer s > 2 and apply the following notation: if A is an s-by-s matrix, then r :— tr(yl) 
and Tfe := tr(^'^) for any integer k > 2. 



Lemma 2.12. Fix a positive integer n > 2 and suppose that A is an s-by-s (componentwise) non-negative 
matrix. Then 

s"-V„ > r". (2.10) 

Moreover, in the n = 2 case equality holds if and only if ak.k = 0,1.1 for k — 2,3, . . . , s and ai,jaj,i = for all 
i 7^ j, further, ifn>3 and s"~^t„ = r", then ak,k — for k — 2,3, . . . , s. 



Proof. If A and B are non-negative s-by-s matrices, then {AB)i, k > cLk^kbk,k- Applying this recursively, 
we get that tr(^") > X]fc=i ^fc fc- Then the inequality between the n"^ power mean and the arithmetic mean 
shows that 

../ tr(A") ^ j ELi«fc,fc ^ ELi «fc.fc ^ 

Vs~\/ s ~ s s' 

proving the trace inequalities. Now if n > 3 and s"~^t„ = t", then we have equality in the power mean 
inequality, which implies ak,k = for A: = 2, 3, . . . , s. Finally, suppose that n — 2. Then st2 = holds if 
and only if we have equality in the power mean inequality and tr(^^) = X)fe=i '^fc the former holds if 

and only if a-^^k — o-i,i for A: = 2, 3, . . . , s, while the latter holds if and only if ai,jaj,i — for all i ^ j. □ 



Remark 2.13. 

• Jf s = 2, n = 3, fli^i = 02.2 — and 01^2 = 12.1 = 1, then s"^^r„ = r", but 01^202,1 7^ 0. 

• The constant s"~^ is the best possible (as shown by A = I). 

• We can not expect a "converse" trace inequality, since for the matrix A with O's in the diagonal and 
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with all other entries 1, we have tr^(A) = and tr{A^) > 0. 



For the conjectured optimal RK method satisfying R{A, b) — 6, we have equality in {2.10) for n = 2. 



Let us now fix any s > 3. By repeatedly using the formulae for the derivative of the determinant and the 
trace 

(det (det$(-)) •tr($-i(-)*'(-)) and (tr $(•))' = tr($'(-)), 

where $ : K — M'*^^ is a smooth — and for the first formula, invertible — matrix function, the cyclic invariance 
of th e trace together with the p — 2 order conditions, we see that the numerator of the stability function 
( 1.1 ), det(/ — zA + zlb^), can be written in the following form 



P(z) 



l + (l-r)z + i((l-r)2-r2) + ^ 



akZ 



k=3 



with suitable real parameters Uk (there may be further restrictions on the parameters which are ignored 

her e). Le t us introduce gq := 1, ai := 1 — t, a2 := — T2). Then, due to A > 0, we have t > and T2 > 0, 

I ... 2 
so (2.10) with n — 2 implies —T2 < — thus 02 

derived a non-linear condition 



|(af — T2) < 5(01 — but r = 1 — ai, hence we have 



> a2. 



(2.11) 



2.5.2. Uniqueness results for polynomials closely related to the class Iis/s.2 In this section, the main 
lemma is Lemma |2.16[ although Lemma |2.14| can be of independent interest, giving lower and upper bounds 
on the coefhcients of a general polynomial P with P(0) = 1 that is absolutely monotonic at a point —r < 0, 
together with a uniqueness result. 



Lemma 2.14. Let s denote a fixed positive integer and set P{z) := 1 + X)ri=i '^n-^" with some real 
coefficients a„. Fix any r > and suppose that P'^^\—r) > for all fc = 0, 1, . . . , s. Then for each index 



1 < n < s we have < 

p{z)^{i + fY. 



< 



— Moreover, if there is at least one n {1 <n < s) with an — 



then 



Proof. By assumption, with 7^ := 



P''')(-r) 



r'^ we have ^ for A: = 0, 1, . . . , s. Taylor expansion, the 



binomial theorem and interchanging the order of summations show that 

k=0 



n=0 fe= 




Now fix 1 < n < s. By equating the coefhcients of z" we get a„ = ^ J2k=n («) 
of the constant terms 



further, from the equality 



k=0 



Ik = 1- 



(2.12) 



Non- negativity of the 7fc coefficients implies a„ > and ( 2.12) sh ows < 7/0 < 1 for all < fc < s. Finally 
we seek the maximum of X]fc=n7fe(^) knowing < 7^ and (2.12). Suppose that 7^ > for some m with 
< m < s — 1. Then we define 7^ := 7^ + 7m, 7^ '■— and 7^ := 74, for s ^ k ^ m. Clearly, < 7^ for each 
< A; < s and X^fc^oT* = 1' ti^t (by using the convention that (™) = if to < n) 



k—n 



1 - E - (^) 7™ + f " ) (-7™) = 7« 

k—7i 



and this last expression is strictly positive, because 1 < n < s is fixed, s > m and 7^, > 0. We can therefore 
conclude that (for fixed s > 1 and 1 < n < s) the value of X]fe=n 7fc(n) O'^®^ all < 7fc (fc = 0, 1, . . . , s) 
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and under condition (2.12) is maximal if and only if 70 = 71 = ... = 7s-i = and 7s = 1. This property 
establishes the upper bound on a„ and the uniqueness part as well. □ 



Next we present an interesting lemma about representing a certain linear combination of the three lowest 
order coefficients of a general polynomial of degree at most s in terms of another linear combination of its 
derivatives evaluated at —2s. 



Lemma 2.15. Fix any integer s>2 and choose an arbitrary P{z) := 1 + X]fc=i '^fe^'^ (with complex 
coefficients). Then 



E 

fe=0 



k\s{s~l) 



s - 1 



Proof. Direct comparison of the am coefficients on both sides by using some simple binomial identities 
(that can be obtained, for example, after (successively) differentiating the binomial expansion of x 1— >■ (1 + 

xY). □ 



The preceding two lemmas are used in proving the following uniqueness result in a class of polynomials 
of degree at most s (unfortunately only in a restricted range s < 8) which are absolutely monotonic at —2s 



and whose lowest order coefficients satisfy a certain non-linear relation. The origin of condition (2.13) is of 



course inequality (2.11) in a slightly different context. 



Lemma 2.16. Let us pick an integer s with 3 < s < 8, and set P{z) := 1 + X]fe=i o-kZ^ with some real 
coefficients a^. Suppose P*^'^)(— 2s) > for all k = 0,1, . . . , s, further, that 



2 , . 

2V"^ ^l-"^- 



ThenP{z) = {l + f^y 



(2.13) 



Proof By Lemma 2.14 with r = 2s, we have a„ > for all n = 1, 2, . . . , s. On the other hand, since now 
p('=)(_2s) > by assumption, by applying Lemma 2.15 we get that 



8s 

1 - 4ai H a2 > 0. 

s — 1 

From this inequality we estimate 02 as 02 > ^^(4ai — 1), so ^ {af 
be factorized as 

(2ai-l) (2ai(s-l)+5-s 



>0, 



- %f(4"i ~ which can 



showing (by using s > 1 only) that either 



a,<l and «i < 



(2.14) 



ai>l and «i > ^^3^- 



(2.15) 



s-5 



We now show that (|2.14| can not occur for 3 < s < 8: indeed, for any 1 < s < 9, we have 2{s-i) ^ 2 ^ 7^ 



< 



i, therefore ai < j. But then (2.13) says that 



- ^ - 2 V s J 2s 

so a^(s — 1) -f 2ai — 1 > 0, implying either ai < < or ai > > ^+^79 ~ i- '^^'^ ^^^^ '^^^^ ruled 

out by ai > 0, and the second one is by oi < 
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Hence the only possibility is (2.15 1 above. But then 5 < ai < 



5 by Lemma 



by the same lemma — we have P{z) — (l + ^) 



2.14 



and so — again 

□ 



Remark 2.17. Notice that ^ contains many uniqueness results for polynomials with maximal radius of 
absolute monotonicity, but polynomials in that article should approximate the exponential function near the 



origin to at least Erst order, i.e., -P(O) = P'{0) ~ 1 should hold. In Lemma 2.14 only P{0) — 1 is required, 
whereas in Lemma 2.16 we assume P(0) = 1 together with the non-linear condition 1^2.13 ). 



2.5.3. The structure ofUg/g p Let us fix 3 < s € N. The non-trivial p values are 2 < p < s + 1, but the 
p = s + 1 case is special, because then Hg/s.p contains only finitely many functions. For 2 < p < s, the general 
form of ip G Tlg/g^p reads as 



(1 - azY 



(2.16) 



with suitable real parameters a,ap+i, . . . ,as, and, of course, with the usual convention that X]m=n(') ~ 0, 
if 71 > A^. Formula (2.2) is recovered if we choose p = 2 in (2.161. 



Remark 2.18. The conjectured optimal ip (i.e., the one with maximal radius of absolute monotoni city) 



in the Tls/s,2 class has o, = whereas in the II^/s 3 class it has a = ^ (^l — ^ fqif^ • in light of Section 
these values have been verified for 3 < s < 4. See |§| formulae (3.1) and (3.2)] also. 



2.3 



Remark 2.19. 



The optimal rational function in the Ils/s,2 class for 1 < s < 4 satisRes the relation 



2.5.4. Possible extension of certain proofs and some questions 



Remark 2.20. In Jl2\ Theorem 4.4], a necessary and sufRcient condition is formulated for a general 
rational function ip satisfying Assumptions ®~@ in our Section \1.1\ to be absolutely monotonic on an 
interval [x,0] C {—B{ip),0]. This condition, among others, requires checking the positivity of the "dominant 
coefRcient" c(ao, /i(ao))i being the numerator in the partial fraction decomposition of ip corresponding to 
the positive real pole ag in DeRnition \1.8\ in our work, and h avin g the highest pole order. It can be seen 
that c(ao,/i(ao)) > is equivalent to our necessary condition (2.3) formulated in the Ils/s.2 class, if "> 0" 



is replaced by "> 0" in {2.3). Moreover, "= 0" in i2.3y holds precisely if the numerator and denominator 



of tp have a common root. We add that we have not used \12\ Theorem 4.4] directly, because the explicit 
construction of its functions F{k,x), L{x) and K{x) (c.f. our Remark 1.13 also) — now possibly depending 
on an additional parameter — would not be straightforward. Instead, we reproduce the "partial fraction 
decomposition + factoring out the dominant term + checking the sign of the remainder in Rnitely many 
cases" idea behind the proof of 12 Theorem 4.4] only after the optimal parameter value a* and the optimal 



radius of absolute monotonicity Rs/s,p (or Rs/s,p) have been conjectured. 



Remark 2.21 (on the explicit computation of tp^'^'^). The sum in {2.4) has a bounded number of terms 
for all k E N, yielding a nice representation of the derivatives of functions in Ils/s,2- For ip e I^s/s,p, a similar 
representation would require a sum with unbounded number of terms as k increases, which precludes the 
extension of our analysis for ■0 e n^/s ^ to the more general class. Fortunately, the optimal ip £ ^[2/2,3 is also 
an element of Il2/2,3- 
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Remark 2.22. In the proof of Theorem 2.7 we wish to show that (under condition a > 0) '\Jj^^\~2s) > 



(k €N) impUes uniqueness. Instead of appealing to formula (2.4) for ■0^'^^ , one could first write up the partial 
fraction decomposition of tp (depending on the parameters a, as, ... , as), then differentiate, and get ip^'^^ in 
a different form. In this form, the "basis elements" 

1, fc + 1, (fc + i)(fc + 2), . . . , n^-ii(fc + m) 

appear naturally in the numerators of the partial fraction decomposition of ^^'^K Some of our observations 
suggest that changing this "rising" basis to the "falling" one 

1, fc - 1, (fc - i)(fc - 2), ... , n^-ii(fc - to) 

could be advantageous: for the optimal ip function, coefficients in the numerators in the falling basis are 
simpler (in fact, for 3 < s < 4, they are 0, 0, . . . , 0, , l_i ) than in the rising basis. We add that the sum 



in {2.4) provides us with a third type of k-basis. We do not know whether this approach could be used to 
characterize the optimal £ ^s/s,2- 



Remark 2.23 (on the point and interval conditions). In our proof of Theorem 2.7 for 3 < s < 4 we show 



that a > 0, (2.6) and (2.7) imply uniqueness. According to our investigations with Mathematica however, 
different formulations of absolute monotonicity of tp on [—2s, 0] can lead to the same uniqueness result. For 
example, if s = 3 and a > 0, then 

ip^''H-(i)>0 {k^ 0,1,2) and & [-6,0] : ip'^^^x) > 

implies uniqueness. Moreover, if s — S, a > and A > 0, then 

i/'('=)(-6) > (fc = 0,l,2) and Vx e [-X,0] : i^^'^'^ (x) > 

implies uniqueness if and only if 0.103209 < A < 10.5836 (the exact values of these algebraic numbers are 
known). We also have the following observations on the "gaps" in the a > parameter region: if s £ {3, 4, 5}, 
a>0 and ip'^^\-2s) > (k = 0,1, . . . , s - I), then 

1 

- 2s 



(recall that for the optimal ip, we conjecture a — as mentioned in Remark 2.18). Moreover, there are 
other gaps simultaneously implied by the above conditions: 

• the interval ^ < a < I, if s = 3, 

• the interval g < a < ^ , if s = 4, and 

• the interval < a < 0.1835 (with the upper bound being the root of a cubic polynomial), if s — 5, are 
also excluded. 

Remark 2.24 (on the Mathematica implementations). This work could not have been completed without 
Mathematica, and especially, without one of its key commands in polynomial algebra. Reduce. Both in the 
conjecture and proof phase, Mathematica's abilities to manipulate high-order root objects symbolically, 
numerically or graphically played a crucial role. Throughout our investigations, we tried to apply symbolic 
methods to the fullest extent possible. At the same time, each pivotal formula was tested numerically with 
several parameter values (some by using different Mathematica algorithms). For typical checks we have 
used 25 digits of precision, but when huge coefhcients were involved, we switched to 100 or 1000 digits 
of precision. Results of the symbolic computations and those given by the numerical algorithms were also 
compared graphically, and found to be in perfect agreement. In fact, as a by-product of our investigations 
we have discovered and reported a few minor Mathematica bugs, occurring with very low probability, in 
contexts when certain algebraic relations between different parameter values are satisfied. 

Let us close this section by posing two natural questions. 



Question. What properties — e.g., what geometrical configuration of the roots and poles — of the rational 
function ijj determine whether R{'4') = B{i(j) or R{ip) < B{tp)7 If R{ip) < B{ip), then which roots of which 
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derivatives of i/' will limit the value of R{ip) in the sense of Theorem [lTi]? 

Question (c.f. the uniqueness of the optimal polynomials in |9|). Let us consider a given {s,p) pair 
(s > l,p > 2). Is there always a unique rational function ip — ^ ^ ^s/s,p (or V' = ^ ^ ns/^ p) with P(0) ~ 
Q(0) 1 such that = Rs/s,p (or i?(-0) = Rs/s,p)'^ 

2.6. Notation used in the proof sections 

For 3 < n E N,n > m E N+ and aj G M (j = 0, 1, . . . , n) the m"^ root of the polynomial equation a„2;" + 
a„_iz"~"'^ + . . . + fliz + flo = will be denoted by 

root„j(a„, a„_i, . . . ,ai,ao). 

We use the following ordering for the roots: real roots are listed first in increasing order, then non-real 
complex roots follow as their real part increases, finally, complex roots with equal real part are sorted 
according to increasing imaginary part (in this context, multiple roots will not be encountered). So for 
example, the value of root^(l, 0, 0, 0, 0, 0, -1) for m = 1, 2, . . . , 6 is -1,1, ~^+'^ , and iif^, 

respectively, corresponding to the equation — 1 = 0. We will give numerical approximations to these roots 
as well, typically to 6 digits, if they occur only in auxiliary computations, but more digits will be displayed 
for significant constants. We will use the « symbol without rounding, i.e., in the sense that all shown digits 
are correct. However, we underline that these approximations are only given for the sake of the reader's 
convenience, and the proofs are not based upon them. 

In our arguments, several exact algebraic numbers (typically denoted by Qspn), polynomials (denoted by 
Pspn) or other auxiliary functions (fspn) will appear. Here the subscripts s and p correspond to the ones in 
n^/s.p or lls/s,p, whereas the positive integer n serves as a counter within the actual section and is increased 
sequentially. Superscript * will often be used to denote optimal values within the family (e.g., a*p being the 
parameter value corresponding to the rational function with maximal radius of absolute monotonicity in the 
class Hg/^ p or H^/^ p). If tf; depends on a parameter, say, on a, then the fc**^ derivative (/c G N) of z i-> ipa{z) is 
denoted by ipi''\ Moreover, we will often use the simpler forms R{a) and B{a) instead of R{i^a) and Biipa), 
based on the bijection a ^a- 

Remark 2.25. We remark that the function a i— > R{a) E [0, +oo] can be discontinuous. We will also see 
examples when the function a i— )• B{a) E [0, +oo] is non-differentiable, or convex on an interval and concave 
on another one. 



3. Determination of R2/ 2. 3 



The functions in H2/2,3 have the form 



(| + f)z^ + (a+l)z + l 



with a e M. For any a € M, the fraction can not be simplified, since the resultant of the numerator and the 
denominator w.r.t. z \s ^Q. This implies that '0a(-) does not have removable singularities, only pole(s). 
This last statement remains valid \i'4>a{ ) is replaced by its A;'^ derivative {k > 1), since differentiation cannot 
introduce new removable or pole singularities, nor annihilate an existing pole. According to Theorem |1.10[ 
R{ipa) > implies that ipa has a positive real pole. This last condition holds if and only if a < —1 — w 
-1.57735 or a > -1 + 4? ~ -0.42264. 

Figure [ij suggests that the maximal radius of absolute monotonicity within this class is 1 + -s/S, 
corresponding to parameter value = + (when the rational function has a pole of maximal possible 
order, i.e., 2). 

For a < —1 — it is easily seen that ipa has only positive real poles. Furthermore, we have ^a(— 2) = 
7Y2^ > and ifia (— |) = — j^^^jf < 0, so by the intermediate value theorem ipa{-) has a root in (— |, —2) 



implying by Theorem 1.11 with i = that for these a values i?(V'a) < | < 1 + V^. 
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Figure 1. The figure shows roots of the derivatives ip\t\-) as a function of a for < fc < 5 for the family 
of functions in the 112/2,3 class. The parameter region between the vertical dashed lines 

^z.e., — 1— ^<a<— 1 + should be ignored, since here R{ipa) = 0. The black dot indicates the 
optimal a parameter value and (the negative of) the optimal radius of absolute monotonicity within this 

class. 



On the other hand, if a > — 1 + > — 1 + then has a pole in (—1 — -\/3, 0), meaning that R{tpa) < 
1 + -s/S here, since ipa is not defined on the whole (—1 — a/3, 0]. Now we are going to exclude the parameter 
■-l + 4.,-l + ^].Ifa=-l 

1 + ^,-1 



region a e 
Let a E 



1 

73 
1 

73' 

and then — by Theorem 



^ j be fixed. Then V'al-l - VS) < < V'a(O), so V-a has a root in (-1 - VS, 0) 



1.11 



Finally, we see that if a 



with £ 

"^23 = " 



+ then R{ipa) < 2 < 1 + ^3, because V'a(-2) < 0. 
hen ^;(-l - a/3) < < ^;(0), 

1 — inequality R{ipa) < 1 + \/3 holds, 
en 

3) ((\/3-l) + 2V3 x + 6) 



/3 + 3-a;) 

from which one can prove recursively for fc > 1 that 

, (fe) ^ _ 6fc! ((12fc - 3) + V3{7k - 2) + (3 + 2a/3) 



k+2 



With the help of these formulae, one verifies that V'!' i^) — fo^' all fc > and —1 — \/3 < x < 0, proving 
our claim. 



4. Determination of Rs/s^p for 5 < p < 6 and a lower bound when p = 2 



4.1. Determination of R^/^fi 



The set 113/3 g consists of the unique member ■036 (^) = — ^"2 — ■ This Fade approximation of the 
exponential function has radius of absolute monotonicity 

R{^36) = rooti(l,-6,0,40) = -2+ s/i ^1^3-75+ \/ 3 + Vb^ w 2.207606. 
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To prove it, we first determine i3(V'36) (see Definition 1.8). Let ao := rooti(l, —12, 60, —120) « 4.64437 and 
a\ := root2(l, —12, 60, —120) « 3.67781 — i ■ 3.50876 denote the real and one of the complex roots of the 
denominator of V'se- Then x ^ x* —B{'ijj^Q) < is the unique real solution of the equation 



|ao - x\ 



\ai 



From this we see that i3(V'36) = rooti(l, — 6, 0, 40). From Theorem 1.9 we know that R{ip26) < B{'ih^ 



In order to show that R{ip3e) = B{tp3Q), we verify that ipsei^*) ^ for any k> 0. (Then Theorem 



1.12 



guarantees absolute monotonicity on [—B(^3q),0], since the denominator of i/j^q does not 



with X = 

vanish on, say, (— oo,0].) The partial fraction decomposition of ■036 is ipaeix) = — 
with Co := rooti(l, -24,-1200, -40000) « 57.2025 and a := root3(l, -24, -1200, 
20.5831. But i^seix*) = rooti(l, -15, 2127, -233) > 0, and for fc > 1 we have 



(fc)/ 
36 \ 



kl{aQ 



N-fe-l 



Co + ci 



ai 



k+l 



ci , 

"1 



ao-x 

-40000) ; 



fc+lN 



+ 



ai—x ai—x 

-16.6012 + ^ 



Since 







ao-x* 


Ctl—X* 




CKi—X* 



1 by construction and 52.8874 w |ci| + |ci| < Co, the positivity of i^^'^ix*) follows 



4.2. Determination of R 



3/3,5 

The set 113/3 5 can be written in terms of one real parameter a G M as 

^ (a+g^)r^ + (6a+|^)z^ + (l2a+|)z + l 

The numerator and denominator do not have a common root for any a, because their resultant w.r.t. z is 
8640 7^ ^- 

Theorem 1.10 tells us that if R{tpa) > 0, then ipa has at least one positive real pole. We claim that the 
+ — 6a) x'^ + (l2a — |) x + 1 = implies a < 0. Indeed, after a rearrangement 



statement 3 x > 
we get 



8x - 20 



20x (x2 - 6a; + 12) 



< 0, 



(4.1) 



since the numerator is negative and the denominator is positive for all a; > 0. 

Next we prove via the intermediate value theorem that for a < — the 10**^ derivative of '0a(-) has a root 
between —2 and 0: indeed, for a < 



' 50 



v^rn-2)<o<v^r^(o). 

The formal proof of this last statement is completely analogous to that of the corresponding statement in the 
114/4,7 class in the opening paragraphs of Section [5[ hence we omit the details here. The only thing we have 
to verify is that the intermediate value theorem is applicable in [—2,0]. On one hand, the poles of ip^^\-) 
coincide with those of V'a(') (only their order can be different). On the other hand, for any fixed a < 0, the 
function V'a(') does not have any poles in [—2,0] (nor in (— oo,0]), because for a: £ R \ {0}, inequality (4.1) 
implies a; > 0. This means that for a < — ^ we have R{ipa) < 2 due to Theorem 1 1 ■ 1 1| with I — 10. 

So in order to get R{4'a) > 2, it is enough to consider a e (— g^, O). As it turns out, the maximal radius 
of absolute monotonicity for these rational functions will be determined by the maximum of the function 
a I-)- B{a) 



B^ipa) according to Definition 1.8 and Theorem |1.9[ 
First we give the real and imaginary parts of the three poles (even for all a < 0) , but due to symmetry, it 
is enough to take into account only the closed upper half-plane. If a < 0, then the real pole is located at 

^351 (a) = rooti(20a, 1 - 120a, 240a - 8, 20). 

As for the other pole, first we set 

£,352 = rooti (864000, 43200, 360, 1) = - ^ (2 -f -I- v^j w -0.0403944. 
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Then we consider the denominator of ipa, substitute z = x + iy, separate the real and imaginary parts and 
solve the resulting system: if £1352 ^ a < 0, then the upper complex pole is found at 




14400a2 + 1200a + 4 - (240a + l)e353(a) 
6(M2^^^^^(^^)y^l 



ith 



g353(a) = rootl(1600a^ -QGOOa^ + 80a, 19200a2 - 400a + 1, -14400a2 + 400a - 4), 



while for a = 51352, the upper complex pole is located at 4-v^ + i (2^/2 + ^/4 - 2) sa 2.4126 + i ■ 2.10724. 
A case separation is necessary in the above formula to avoid a 0/0 under the square root: geometrically this 
a = 5)352 value corresponds to the case when the three poles have the same real part. 




Figure 2. By substituting z = x + iy (x,y £R) into the fraction in {4.1) and separating real and 
imaginary parts, two M surfaces (Src and Si^) are obtained. The contour Sjm = yields the equation 

y{x'^ — 16x^ + 2x^y'^+ 96x^ — 16xy^ — 2A0x + y^ + 16y^ + 240) = depicted as red curves (with three 
connected components). They describe the locus of the poles of 4'a for a € M (including the limiting case 
a — ±00) on the complex plane in the 113/3 5 class. For a particular, fixed a € M, the three poles of ipa are 
found when the red curves intersect the contour curves Sre = a (some of which are shown as blue curves). 
The poles in the a = ±00 limiting case are located near the three singular points (squeezed in between the 

white regions), where Src and Sim are undefined. 



Now with the real and imaginary parts separated, it is easy to find, for any a < the unique point on the 
real axis being equidistant from the real pole and the upper complex pole. We check that for a < 51352 the 
(positive) real pole is strictly smaller than the real part of the (upper) complex pole, hence — by applying 



the notations of Definition 1.8 with ao = ao(a) being the positive real pole — the equidistant point on the 
real axis is strictly positive, meaning that £ J(ao(a)) and B{a) = — inf J(Q;o(a)) = +cxi here. We easily see 
that i?(5'352) = +00 also, but for 51352 < a < 0, B{a) will be finite. If 51352 < a < 0, then the point on the real 
axis equidistant from the real pole and the complex pole is found at -^^"l^j with /351(a) :— 

glr,,{a) (60a5?353(a) - 120a + 1) - SOaglr^ " (2 - 2i0a) gl^-^ia) - (240a - 8)5^353(0) - 20 

and 

/352(a) := 2 (5»35i(a) - 5'353(a)) (60a5»353(a) - 120a + 1) . 



/351(a) 
/352(a) 

5,354 = rooti (13824000, 345600, 2880,-152,-1) « -0.00625485, 



Now Reduce establishes quickly that has a unique zero in 51352 < a < 0: by defining 
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<o, ^ 

(5352, ^354) J 352 



_ Q /35I 

a=Q35i /352 



> 0. 

(e354,0) 



it turns out that 

/351 

/352 

(However, our "manual" computations to reproduce this sign property without Reduce occupy several pages 
in the corresponding Mathematica notebook when written out fully, so we only sketch the main steps. From 
the equation /351(a) = we simply express £'351(0) as £'351(0) = _[_ V2-\/40ag353(°) - "^ where we have displayed 

-^60ae353(a)-... 

(and will display) only some (typical) "highest order" terms to indicate some structure. This ± expression 
is then substituted back in place of x into the defining polynomial equation 20ax^ + (1 — 120a)x^ + (240a — 
8)x + 20 = of the root object ^351 (a). After clearing the denominators, we arrive at a relation of the form 



T(40a(120a - 1)91^ - . . .) v/60ag353(a) 
732 • {200a'gl,,ia) + . . AOagl,,{a) - 



with the square roots being the same as above in the ± expression. A squaring and rearranging cancel the =F 
symbol and the square roots, and yield a polynomial in a and in £3353(0) of the form — 51200000a^£'353(a) + 
... + 80 = 0, with 39 terms omitted in between. Now we express £'353(0) from the defining polynomial 
equation of this root object ^353 (a) as £1353(0) — 

— r ((9600a2 - 80a) £'353(0) + (-192000^ + 400a - l) ^353(0) + 144000^ - 400a + 4) . 
loOOa^ 

This formula is applied recursively within the above equation — 51200000a^£'353(a) + . . . + 80 = to eliminate 
all the powers £1353 (a) with fc = 9, 8, . . . , 3 and to have 

-f353(a) + ^354 (0)^353 (a) + -P355(a)£'353(a) = 0, 

where ^35^ (a) (£ — 3, 4, 5) are suitable polynomials in a of degree 7 with integer coefficients having at most 
14 digits. The above quadratic polynomial is now solved to get 



, , _ -P354(a) ± yPlija) - 4^355(^)^353(0) 
£'353 (.a j - — — ^ , 

2^355(0) 

which information is then substituted back into the very defining polynomial equation of the root object 
^353 (a), and we can derive the equation P356(a) ± -P357(a)\/-P358(a) = 0. (Here P356(a) is a polynomial in a 
of degree 18 with integer coefficients having at most 35 digits. The degree of the integer polynomial ^357(0) 
is 10, while that of ^353(0) is 14.) A final squaring and rearrangement allow a factorization of the form 

= a (864000a^ + 43200a2 + 360a + l) (l3824000a^ + 3456000^ + 28800^ - 152a - l) x 
(I7280000a^ + 1728000a^ + 576000^ + 440a + l) x 

(248832000000a^ + 26956800000a^ + 6393600000^^ + 14112000a^ + 348000^ - 310a - l) . 

Now we see that only 3 out of the 18 roots of the above polynomial are located in the (real) interval (£1352, 0), 
and a check (similar in spirit to the manipulations so far) yields that 2 out of these 3 are not solutions to the 
equation /353(a) = 0, just the third one. (Of course, we have checked during all the computations that the 
quantities appearing under the square roots are non-negative, further, that the possible vanishing leading 
coefficients of the intermediate linear or quadratic equations do not yield extra solutions to /353(a) = 0.) 
Summarizing the above, we have proved independently of Reduce that the equation /351(a) = has a unique 
solution £1354 in (^352,0), and a continuity argument yields also the signs of the fraction on the 

subintervals.) 

Now let us proceed further. Knowing the sign of implies that on a G (£'352,0) \ {£'354} relation 

e /(Q!o(a)) can hold only for a G (£1352, £'354), so 

B{a) = 0, if a e [£1354, 0), 

and 

B{a) = -inf /(ao(a)) = - {^^^^°| , if a G (^352, £'354)- 

/352(a) 
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Figure 3. In the n3/3_5 class the function a i-^ —B{a) is shown as the thick red curve (with = —oo to 
the left of the long, vertical dashed line). The small dashed rectangle will contain the optimal value of —B. 

The other curves are the roots of the derivatives ^^\-) for fc = 0, 1, 2, 8, 9, 10, 11, 12, 13. In this plot 
window, these are the only roots of 4''a\') for k < 13. We remark that altogether 17947 digits are needed 
just to write down the integer polynomials appearing in the numerators of 'ipa'\') for these k values. The 
structures analogous to pitchfork bifurcations in this figure nicely illustrate RoUe's theorem, further, the 
fact that if a smooth function has a root of multiplicity m at a point, then its derivative also has a root 

there (with multiplicity m — 1). 



Having completely described the function B, in the next step we wiU determine its maximal value on 
a e (~ O) by differentiation. We check by computing the discriminant of the denominator of \pa that all 
roots of the denominator are simple for all a < 0. Moreover, this denominator is smooth in x and a, so 
the implicit function theorem yields that the functions a £1351(0) and a 1— >■ £?353(a) are differentiable with 
derivatives 



^ 20e3^i (g) - 120g2^i (g) + 240g35i (a) 

^^^^^ ' 60gg^5i(g) + 2(l-120g)£.35i(a) + 240a- 

£•353 (a) 



8 

and 

3200g£»i53(a) + (-19200a + 80)ei53(g) + (38400g - 400)g353(fl) - 28800g + 400 
4800a2g2g3(g) + 2 (-9600g2 + 80g) 5353(0) + 19200g2 - 400a + 1 ' 

Using also these formulae we determine that there are exactly two a values in (— O) such that B'{a) — 0, 
one is a local minimum and the other one is a local maximum of B. But since we are maximizing B, the 
local minimum is ignored. The local maximum of B is located at 

«35 := 

roots (13436928000000, 1492992000000,68428800000, 1427328000, 13867200,61200, 101) w 

-0.009257142292762484937472363 

with value 

B{al^) = root2(l, 12, 36, -76, -360, 0, 900) « 2.301322934003485801187482. 

Note that Figure [3] depicts —B so that it can be compared to the derivatives of V'o(')- We remark that 
the local maximum of — q-j is described by the same root object as agg but with roots replaced by 
rooti. Reduce was able to solve the equation B'{a) = within 2 seconds. When we solved the same equation 
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"manually" (i.e., in a Reduce-independent way) with all the details written out fully, the corresponding 
Mathematica notebook saved as a PDF file consisted of approx. 280 pages, and the intermediate expressions 
could be described by approximately 700 000 characters. The final integer polynomial with both root objects 
£•351(0) and ^353(0) eliminated has degree 162 and leading coefficient approx. 4 • 10'^^^. In view of the above, 
it is very surprising that the algebraic number 035 is so simple: the reason is that the final polynomial of 
degree 162 can be written as — 64000a** P35g(a)P35iQ(a)P3^|2i('^)^35i2('^)-P35i3(o) '^ifh the polynomials P^^e 
(£ = 9, 10, . . . , 13) having degrees 10, 26, 3, 4 and 6, respectively, further, ^3513(035) = 0. 

As a final step, we show that the interval [—B (ogg) , 0] does not contain any roots of the derivatives ip^? (•) 
(A; = 0, 1, . . .). In view of Theorem 1.12 with x — —B (ftgg), it is enough to show that for all < A: € N we 
have 

V'i2(-^K5))>o, 

which we check by partial fraction decomposition (also showing that there are no poles of ipa* in [~B (035) , 0] , 
applies). The function ijja* evaluated at x admits the following decomposition: 



so Theorem 



1.12 



Co 



Cl 



Cl 



— X oil — X ai — X 

with some c, Cq G M, ao > and Ci, ai G C \ M. These constants are algebraic numbers of degree not exceeding 
12, and absolute value of the coefficients of their defining integer polynomials less that 10^^. Here we give 
only their numerical approximations as 

c « -0.800411, Co « 46.829419, ao « 4.449434, 



Cl w -13.750016 - i ■ 16.640148, ai w 3.475900 + i ■ 3.492336. 
With X = X* := —B (035) < we check that 

0.09999997 « ?/'a*(a;*) > 0. 



On the other hand, if fc > 1, then ip^J (x*) takes the form 



fc!(ao — X*) 



*\-fc-i 



Co + Cl 



/ ao - X* 
\ai — X* 



k + l 



Cl 



ai ~ X* 



fc+i^ 



Now we make use of the facts that ao 



> and 



Co > and 43.172096 w |ci| + |cr| < cq. This proves that V^-^ (x*) > also for fc > 1. 

Summarizing the above, the optimal radius of absolute monotonicity within the 113/3 5 class is B (ajg) 
2.301322. 



= 1 by construction, further, that 



4.3. A counterexample to the conjecture that Rs/s,2 = 2s 

A counterexample to Conj ecture 1 1 . 1 1 for the case m — n = 3 was found via extensive numerical search using 
MATLAB and the Symbolic Toolbox employing the fminsearch function (Nelder-Mead simplex method); 
it originally appeared in [7j . For the sake of our presentation, let us denote this counterexample with rational 
coefficients by '032 Now we are going to present a simpler counterexample — in the sense that the rational 
coefficients have much smaller numerators and denominators — with slightly improved radius of absolute 
monotonicity. ^ 

By suitably perturbing the rational coefficients into nearby simpler ones and embedding into a one- 
parameter family of rational functions in a way that the order conditions are satisfied within the family, say, 
as 



2289 2 
34970 



269 



Z + 1 



327 



150 
269 



Z + 1 



tWe note that there is a typo in the last digit of one of its coefficients in |7j, instead of 7969150767159903, one should have 
7969150767159904. 
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with c £ M, we can optimize the radius of absolute monotonicity w.r.t. c using Mathematica. It turns out that 
the optimal parameter c* (i.e., the one that yields the maximal R within this chosen class) is an algebraic 
number of degree 5. By replacing c* with a nearby simple rational number we get, for example, 



1246 ,3 



384649 



34970 



269 



Z + 1 



327-* 



_8_72 
65"^ 



150 , 
269^ 



This function has 



i?(Vj32) = rooti (43572620, -880461561, 5950520030, -13451175530) 



6.778307398562974637718719 > 6. 

We remark that 6.77823595 « i?(V'32) < i?(V'32) < i?(V'c-) ~ 6.77830907. 

Let us briefly give some details that can be used to verify the above value of R{ip32)- The function tp32 at 
X has the following partial fraction decomposition 



Co 



ci 



ci 



with c 



769298 



Co w 24.8122, ao 



"0 - 
3.71417, 



X 



ci 



ai — X ai 

s -8.39838 - i ■ 9.53528 and ai 



3.17368 + ^-3.45514. 

(These are all algebraic numbers of degree 3, and with maximal absolute value of the coefficients in their 
defining integer polynomials appro ximat ely 1.046 • 10^^.) We see that "032 has a unique real root at x* := 
-R{ip32) ~ -6.778307, so Theorem \lu\ with x = x* and £ = applies. We claim that ^"32^ (a;*) > holds 
for all A: > 1 as well. If A; > 1, then 



^-32' i^" 



)^k\iao-x*) 



-fe-i 



Co + ci 



ao — x 
ai — X* 



k+l 



Ci 



x^ 



fc+lN 



Now since ao 



> and 



Qi— CKi— a 

integers as coefficients), we see that the sufficient condition 



0.995991 (an algebraic number of degree 18, with 92-digit 



2|ci| 



ao 



ai 



fe+i 



< Co 



for the non-negativity of '032'' (^*) holds for fc > 5. Finally, we directly check that 032'' (a;*) > is also valid for 
1 <k < A. (We remark that now B{ips2) ~ 7.59982 is an algebraic number of degree 3, so for this particular 
function Riip-32) < Absolute monotonicity in the whole [a;*,0] interval is guaranteed by Theorem 

1.12 by taking into account that there are no poles of "032 hi [x* , 0]. 

Of course, it is to be emphasized that the maximal radius of absolute monotonicity in the whole 113/3 2 
class is still unknown: this class of rational functions can be described by 4 parameters. It is a major open 
challenge to find the maximal R within this 4-parameter family. 



5. Determination of -R4/4.7 
The set 114/4 7 can be described by one real parameter a e M as follows 



^{7a + 4:)z'^ + 2To(21a + 13)^3 + (-i(7a + 2) + a + i) + (a + l)z + 1 
-gioC^" + 3)^4 + 2To(21« + 8)^3 - j^(7a + 2)z^ + az + 1 ' 



The numerator and denominator do not have common factors for any a, since their resultant w.r. t. z i s 
870912000 7^ order to have R{ipa) > 0, ipa needs to have at least one positive real pole (Theorem 1.10 1. 
To obtain some preliminary information on the location of the real poles of a; 1—^ ipai^), we apply the same 
trick as at the beginning of Section 4.2 This time the denominator of ipaix) (with x € M) vanishes if and 
only if 

_ -3.t4 + 32x^ - 120a;2 + 840 
" ^ 7x (a;3 - 12x2 ^ gOx - 120) ' 
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Then, by analyzing the range of this rational function on the right-hand side (under the restriction a: > 0) , 
we prove that ipa has at least one positive real pole precisely if 

a < 5471 ~ -0.843194 or a > 5472 ~ -0.471357, 

so in the rest of this section we can consider only the above parameter set. We remark that the exact values 
of these and the following algebraic numbers are listed in the Appendix. Furthermore, with this technique 
we also establish that for any a < the denominator of ^pa{■) does not vanish in [—1,0]. 

Now we show that in the left unbounded component, that is for a S (—00, 5471], we have R{a) = R{'4>a) l£ 1- 
To this end, let us consider the 9*^ derivative of x 1— ipa{x), which can be written as 

P47(a;, a) ■ 1451520 ((7a + 3)a;'^ + (-84a - 32)x^ + (420a + I20)x^ - 840ax - 840)"^° , 

where ^47(0;, a) is a polynomial in x (of degree 30) and a (of degree 10), and 5917 digits are needed to write 
down all of its integer coefficients. Since, as we have just seen in the previous paragraph, the (—10)**^ power is 
never singular for a < 5471 < and a; € [—1,0], we can apply the intermediate value theorem in the interval 
[—1,0]. On one hand, 

P47(0,a) = 227219566849228800000000 - 433782809439436800000000a2 < 

for ag (—00,5471]. On the other hand P47(— 1, a) = ^471(0), a polynomial of degree 10 with integer 
coefficients. But ^471^ is a positive constant, and for fc = 9, 8, . . . , one computes that sgn ^^471 (—5)^ = 

(— 1)*^', showing inductively that for any a e (—00, — |] and k — 9,8, ... ,0 we have sgn (^Pjl^l (a)^ = (— 1)*"' 
(notice that we have used only rational arithmetic). Taking into account (—00,—^] D (—00,5471], we get 
that ^47i|(-_o2 g^^^] > 0. The intermediate value theorem now guarantees that for each a e (—00,5471], the 

function ip(P {■) has a root in [-1,0]. Due to Theorem [1T1I with £ ^ 9, we have R{il)a) ^ 1 for a G (—00, 5471]. 
Next we focus our attention on the right unbounded component a £ [5472, +00). In this region we are going 



to explicitly compute the function a ^ B{a) = B{ipa) (according to Definition 1.8 1. 

As a first step, Reduce was able to give a complete description of the locus of the poles in the complex 
plane with real and imaginary parts separated as a is varied. If a = 5472, we have a real pole of order 2 
at 5473 « 7.64527. For 5472 < a < — |, it bifurcates into two real poles (of order 1) located at 54711(a) e 
(5474, 5473), and 54712(a) S (5473, +00), with 5474 w 5.64849. If a = — |, then there is a unique real pole of 
order 1 at 5474. Finally, for a > — | the smaller (negative) real pole is located at 54711(a) < 0, while the 
larger (positive) pole is found at 54712(a). 

As for the complex poles, due to symmetry, it is enough to describe only the one with imaginary part > 0. 
For a S [5472, — f ] , the complex pole in the upper half-plane has real part 5477(a) and imaginary part 



/ (7a + 3)5^77(0) - (63a + 24)5177(0) + (210a + 60)5477(0) - 210a 
y (7a-t- 3)5477(a) - 21a - 8 ' 

If a > — I , then the expressions for the real and the imaginary parts of the upper complex pole are obtained 
as above, but with changing rooti in the definition of 5477(a) (see the Appendix) to root2. 

Figure |4] depicts the locus of the poles of -00 as a is varied. For the sake of completeness, we have also 
included the poles corresponding to the interval a G (—00, 5472). 

Now we combine the above formulae on the locus of the poles, real and imaginary parts separated, 
to determine B{a) for each a G [5472, +00). With the help of Reduce (and by applying some geometric 
reformulations to be able to obtain the results in a reasonable amount of computing time, but now omitting 
the details here) we can prove that for —0.471357 « 5472 < a < 5478 « —0.469514 the point on the real axis 
equidistant from the smaller (positive) real pole and the upper complex pole is positive (i.e., lies strictly 



closer to the non-real upper complex pole for a G [5472, 5478)): hence (in the sense of Definition 1.8 1 ^ I{oiq), 
so B{a) = here. The larger real pole — existing for a G (5472, — f ) — clearly cannot infiuence the value of 
B{a) now. For a G [5478, — |] , the equidistant point on the real axis from the smaller positive real pole and 
the upper complex pole is non-positive, and is given by with 

/471(a) := 2(7a + i)gl'j'j{a) - 4(21a + 8)5477(0) + 30(7a -I- 2)5477(0) + 
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6p 




-6 - 



Figure 4. The trajectory of the four poles of ipa in the complex plane in the 114/4 7 class as parameter a 
traverses through the real line from —00 to +00. As a ^ — 3/7~ and passes through this point, the larger 
real pole is repelled to +00, then it enters the real line again from ~oo and approaches to 0^. (In order to 
create the figure, we have used an equidistant grid on a € [—3, 3] and Mathematica's NSolve approximated 
the roots of the resulting 6001 quartic equations appearing as denominators in a few seconds.) Compare 
these curves with the corresponding red curves in Figure [2[ 



(21a + 8 - (7a + 3)£<477(a)) glmia) - 210a 

and 

/472(a) := 2 ((7a + 3)£»477(a) - 21a - 8) (£1477(0) - 6'47ii(a)) ■ 

(Of course, as already noted earlier, the expression "smaller positive real pole" in the previous sentence should 
be interpreted as the "unique real pole", if a = — |. Moreover, the expression "non-positive" can be replaced 
by "zero" if a = g478, and "negative" if a S (^473, ~f ] ■) Finally, we study the interval a > — |. Now one should 
take into account the newly created negative real pole as well. It can be proved that if a = £1479 ~ —0.358565, 
then there is a unique negative real number, rooti(l, —21, 165, —520, 0, 3600, —6000) « —2.40614 such that 
it is equidistant from all the four poles. If a € 9479] i then the negative real pole is still too far on the 
left to have an effect on B{a), so for these a values B{a) is obtained as the absolute value of the point on 
the real axis equidistant from the positive real pole and the upper complex pole. Here B(a) will be given in 
terms of /473(a) and /474(a). Expression /473(a) is defined just as /471(a), but with each rooti occurring in 
the definition of ^?477(a) replaced by root2, further, with £14711(0) replaced by £14712(0). Expression /474(a) is 
obtained from /472(a) via the same two replacement rules. If a > £1479, then the 

(equidistant point on the real axis from the positive real pole and the upper complex pole)i 

is smaller than the 

(equidistant point on the real axis from the positive real pole and the negative real pole) 2, 

hence here B{a) is determined solely by the two real poles. Notice that quantity (...)2 is simply the midpoint 
between the two real poles. For 

-0.358565 w £<479 < a < £'4710 ~ -0.274796, 

the expression (...)2 is non-positive, so its absolute value gives B{a). But if a > £'4710, then (...)2 is strictly 
positive, so ^ ^("o) again, thus B{a) — here. 

We now summarize the above information in one formula for B(a). (For completeness' sake we also provide 
some additional information on B{a) in the interval a € (—00, £1472] already investigated earlier. We remark 
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that a = g475 « —0.850052 corresponds to the geometric configuration when ipa has three poles — one positive 
real and two complex poles — with equal real part, and the 4**^, real pole lies to the right.) We have 



B{a) { 



= +00 
> £'476 

= 
= 

/471(a) 
/472(a) 
/473(a) 



if a e (-00, £1475], 

if a e (£'475, £'471], 

if a e (£'471, £'472), 

if a e [£1472, £'478], 

if a e (£'478,-f] , 

7WS) if ae (-f,£^479] , 

5 |£'47ii(a) + £'4712(0)1 if a e (£'479, £'4710] , 
, if a G (£'4710, +00) 

with £1476 « 21.5907. 

After describing the function B, we find its maximal value on a g [£1472, +00), or, equivalently, on 
[^?472, £"4710]- Mathematical s Maximize was finally able to locate this unique point. It turns out that 
B\, , attains its maximal value at 

I l£?472 ,i?4710j 

047 



-0.4398493860002001824004494, 



where a^j is the smaller real root of a polynomial with integer coefficients and of degree 30. The corresponding 
maximal value is 

B«7) w 2.743911895676330804848228, 
expressible as the larger real root of a polynomial with integer coefficients and again of degree 30. 



Figure 5. The graph of the function a i-^ —B{a) (red, thick curve carrying the black dot) for £(472 < a < 
is depicted together with the roots of the O"^ and 1^' derivatives of ipa(:) in the class. B{a) — also 

between the two vertical dashed lines on the left. The distance between the black dot and the horizontal 

axis is the optimal radius of absolute monotonicity. 



As a last step, we check that no roots of the derivatives of the function ^ai.^i') can enter the interval 
[—B (a|7) ,0], so its radius of absolute monotonicity is as large as it can be (Theorem 1.9), based exclusively 
on the location of the poles: in view of Theorem 1.12 with x = —B (047), it is enough to show that for all 
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< e N we have 

Again, the most convenient form to use is the partial fraction decomposition ofjpaU^x) (this decomposition 



also reveals that 'ipa' has no poles in the interval [—B (047) , 0], hence Theorem 1.12 is applicable). However, 
in this case, due to the high degree polynomial involved in the definition of 047, Mathematica could not find 
the exact partial fraction decomposition of i^) in ^ reasonable amount of time (i.e., expressing the poles 
and the coefficients of the partial fractions as explicit and exact algebraic numbers). In order to overcome 
this difficulty, we have employed complex interval arithmetic with rational endpoints. We started from good 
enough lower and upper rational bounds on a'^^j, i? (a^^) and the 4 poles: we applied validated numerical 
algorithms (such as Isolatinglnterval, with tolerance 10"^*^) in the case of real algebraic numbers, and 
higher precision evaluation in the case of complex roots, which have been previously shown to be the unique 
roots within larger rational rectangles in the complex plane by Isolatinglnterval. Then we also expressed 
the coefficients of the partial fraction decomposition in terms of the poles and a^j symbolically in advance. 
At the end, we were able to give rigorous lower and upper rational bounds on the (absolute value of the) 
coefficients using only rational arithmetic without any difficulty. These computations however produced 
quite lengthy outputs, since the numerators and denominators of some intermediate rational numbers in the 
bounding intervals consisted of integers with more than 220 digits. With the above simple interval technique 
we were able to completely reproduce the numerical partial fraction decomposition of ipai^{x) obtained 
directly in a much simpler manner. 

The partial fraction decomposition of ipai^{x) has the following form: 



aQ — X ai — X a2 ^ X 0.2 — x 

with some c, Cq, ci e M, ag, ai > and C2, 012 G C \ M, further, the following numerical approximations are 
valid (to simplify our presentation, we now omit listing any exact rational bounds): 

c w 11.666934, Co w 202.617318, ~ 5.779490, 
ci « -1084.668490, ai « 47.331517, 
C2 ~ -19.880058 - i ■ 86.974803, 02 « 4.778363 + i ■ 4.007962. 



If fc > 1, then ihi'l' (x) can be written as 

"47 ^ ' 



I ao ~ X 




k\{ao-x) M co+ci 

\ \ai — X 

Now we evaluate the above expression a,t x = x* :— —B {a\-j) < 0. We see by construction that 

— 1, and it can also be proved that < < |. So by taking into account ao > 0, —x* > and 

Co > 0, a sufficient condition for the positivity of vi«'(a^*) {k > 1) is that |ci| • (g)''^^ +2|c2| • 1 < co, or, 
by using the rigorous bounds |ci| < 1085, |c2| < 90 and cq > 202, a weaker sufficient condition is given by 
1085 • (1)''"*'^ + 180 < 202, which is seen to hold if fc > 2. But for fc = and fc = 1 one checks directly (with 
simple rational bounds on a^^ and —B (047), and even without using partial fraction decomposition) that 
^S(-S(al7))>0. 

The above computations confirm that the optimal radius of absolute monotonicity within the 114/4 7 class 
is B(a|7) w 2.743911. 



6. Determination of R2/2.2 and R^/z.p for 2 < p < 4 

As an introduction to the Tls/s,p class, we first list the two optimal rational functions in the TI2/2.P family 
for p = 2 and p = 3. 

In the former case, for the sake of completeness, we also provide two direct proofs of the equality i?2/2,2 = 
4. In the beginning of Section |2.4[ we have already indicated that this equality simply follows from the 
corresponding result in |12) . Our first proof will be similar to the earlier ones — carried out in one-parameter 
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families of rational functions, and not knowing the optimum in advance; while the second proof is based on 
the same idea that we will apply in the Iis/s.2 cases with 3 < s < 4 — having more than one parameter, but 
with an already conjectured optimum. 

Lemma 6.1. R2/2.2 = 4. 



Proof (first proof). Any element of 112/2,2 can be represented as 



\ (2a2 - 4a + 1) z2 + (1 - 2a)z + 1 



(l-az)2 

with a G M. First we perform the same preliminary reduction as we did in Section |2.4[ If a = 0, then the 
rational function reduces to the second degree Taylor polynomial of the exponential function around 0, which 
has i? = 1. So we can assume a 7^ 0. The numerator and denominator have a common root if and only if 



a = ^, in which case i? = 2. If a 7^ i, then due to Theorem 1.10 a > is necessary for i? > 0. So we can 
consider only the parameter values ^ <a^ \. Ifa> 1, then 

^W(O) = -12a2(4a - 3) < < = ^i'^-l), 

mean ing that there is a root in [—1, 0] of the continuous function for these a values, i? < 1 by Theorem 



1.11 



jecause ipa{-l - v^) = 0. If a G (0, 1) \ {|, 1 - ^} 



with £ ^ 4. We will find the optimal R in the interval a G (0, 1) \ {^}, determined by the intersection 
of certain roots of ipa and '(/'a- (Let us add that a plot of the roots of the lowest order derivatives of ipa{') 
in the region < a < 1 would be very similar to Figure 7|)Ifa = l— then the leading coefficient of the 

numerator of ipa vanishes, and in this case 
and a > 4, then one negative root of ipa{') is 

-1 + 2a + v/4a - 1 

> —4 

2a2 - 4a + 1 ' 

while if a G (0, 1) \ {|, 1 — ^} and a < j, then the only root of ^^(•) is > > —4, implying that for 

a ^ J, R < 4. Finally, if a = then ^paiz) = -7— ^4^7 and one directly checks that now i? = 4. □ 

Proof (second proof). We can again consider only functions of the form ^aiz) = — '^'^^i-lz)^^ 2a)2+i 
with a > 0, but now, as opposed to the First proof, we conjecture i?2/2,2 = 4 in advance. We show that 

^a(-4)>0, V;(-4)>0 

and 

VfcGN,fc>2: Vi''H-4)>0 
imply a ~ ^. Indeed, the first two conditions with a > amount to 

IGa^ - 24a + 5 > 0, 12a - 3 > 0, 
that is, to a G {1/4} U [5/4, +00). Then, instead of the third condition, it is sufficient to require the weaker 



inequality (2.3) — expressing the fact that the leading coefffcient of a polynomial in k has to be non-negative. 



if the polynomial is non-negative for all k > 2 — that reads as | — a > 0, implying the unique value a = j. □ 

As for the 112/2.3 class, it has only two members: 

-i(x/3 + l) ^'--^5 + 1 

V'23l(^) = — ,-, ,2 

(l-i(3 + V3).) 

and 



'4'232{z) 



(l-i(3-^/3)z)^ 
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One easily checks that i?(^23i) < i+Vs ^ 5 shown by a root of ^231)- On the other hand, 

R{'^p232) = 1 + a/S, smce ^232 = "^a* ; J-C, the optimal element of the 112/2 3 class, so -R2/2 3 = ^2/2 3 = 1 + 

Vs. 



6.1. Determination of 4 

Any element of 113/3 4 can be represented in the form 

a^z^ + a2z'^ + aiz + 1 
^^'^ ^ [l^^ 

with suitable real parameters a, 01,02 and 03. One directly computes that there are exactly three solutions 
(a, fli, 02, 03) G to the system 

V'(''H0) = 1, fc==0,l,2,3,4 
(we remark that we would get the same three reaJ solutions if we allowed (a, oi, 02, 03) S C*). For all of these 



three solutions, condition a > is automatically satisfied (c.f. Theorem 1.10), and we have for to = 1,2, 3 
that a = rootm(24, —36, 12, —1), ai = root4_TO(8, 12, —12, 1), 02 — root 1 (3„j2_ii„+i2) (64, —48,0, 1) and 03 = 



root] 



(3m2-llm+12) 



]^2)(1536, — 1152, — 24, 1). Let us denote the corresponding rational functions by V'34m ("i 



1, 2, 3). We can check that ^ is a pole of order 3 in each case, so no cancellation between the numerator and 
denominator occurs. We are going to prove that 



max i?(V'34m) = -X* 
m— 1,2,3 



with 



X* rooti(l, 12, 60, 120, -144, -1152, -1536, 1152, 2304,-1536) « 

-3.287278451851993925371346 

and the maximum occurring at to = 1. 

Let us consider the to = 2 case first. Then -0342(0) = rooti(8, -3540, 600, 4625) < 0, so i?(-0342) = by 
definition. 

Next if TO = 3, then V343 (root5(l, 66, 1242, 7008, 1872, -648, -72)) = 0, so Theorem 
that i?(V'343) < -root5(l, 66, 1242, 7008, 1872,-648,-72) « 0.0943315 < 
Finally, if to = 1, then the partial fraction decomposition of -0341 reads as 



1.11 



with £ — 5 says 



0341(2;) = c + ^ 



(ao - xY' 



wherec = rooti(l,9,-9,-9) w -9.82294, ci = root3(l, -288, -5184, 13824) w 304.856, C2 = rooti(l, 3168, -243648, -23500? 
-3243.10, C3 = root3(l, -11808, -684288, 110592) « 11865.6, and ao = root3(l, -12, 36, -24) « 7.75877. 



Then ■034i(a;*) = 0. Since now trivially i?(034i) = +00, Theorem 1.12 with x :— x* proves absolute 



monotonicity of ^341 on [a;*,0] provided that we show the point conditions 

7^f4l(x*)>0 

for all A: > 1 (the theorem is applicable due to the lack of poles in (— oo,0]). Indeed, for such k values we 
have 



4iii^*) = ik+mc^o-x*) 



C3 C2 

2 fc + 2 



(ao - X*) + 



Cl 



-(ao - a; ) 



(/c + l)(fc + 2) 

Clearly, because of ao — > 0, it is enough to verify that the [. . .] expression is non-negative, but since 
limfc_j.+oo[- ■ ■] = ^ > 0, only finitely many k values should be checked. We easily see that [...]> holds if, 
for example, 



ao — X 
k + 2 



|C2| 



-(ao-x*)) <^(3250 
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Now the second inequality above is satisfied for A; > 6, and we determine directly that [• • • ] > holds for each 
k = 1,2, ... ,5 as well, completing the argument that the maximal radius of absolute monotonicity within 



the n 



3/3,4 



class is 



We add that partial fraction decomposition in the m = 
coefficient" C3 is negative, immediately implying R{'ipa2) 



2 case shows that the corresponding "dominant 



(c.f. Remark 2.20) 



6.2. Determination of R3/ 3 ^3 

We will prove that the maximal radius of absolute monotonicity in this 113/3 3 class is 2 
Now the one-parameter family of rational functions takes the form 



/8 « 4.82842. 



ipa{z) 



i {-6a^ + 18a2 - 9a + l) + i (Ga^ - 6a + l) + (1 _ 30)^ + 1 



(l-az)3 

with a e M. If a = 0, then the rational function reduces to the third degree Taylor polynomial of the 
exponential function around 0, which has R — 1. So we can assume a ^ 0. Let us exclude two more exceptional 
parameter values as well. It is easily seen by computing the corresponding resultant that the numerator and 
denominator have a common root if and only if a = g (3 ± -^3) ■ In any of these cases, ipa is either '0231 or 
"0232 from the 112/2,3 class (see the paragraph just before Section 6.1 1, hence R < 1 + \/3 < 2 + a/8- 



Now we can apply Theorem 1.10 saying that a > is necessary to have i? > 0. So in the following we can 
suppose that < a ^ ^ (3 ± VS). 

One can prove by induction that for fc > 1 



(2a - l)(6a - l)a^x^ + {Sa^ik - 2) + 0^(9 - 7k) + a{k - 1)) x+ 



a^(k - 3)(fc - 2) - a{k - 3)(fc - 1) + -(fc - 2)(fc - 1) 



(6.1) 




Figure 6. Tie n3/3_3 class. Each curve (corresponding to different k values between 0, 1, ... ,36 j shows a 
root of 'ipa^\') as the parameter a is varied. The red curve (k = 0) is obtained as a solution of a parametric 
cubic polynomial, whereas the other curves are derived from quadratic ones. 
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Let us consider first the k — 5 special case of the above exphcit formula with a > ^ . From this we see that 
one of the roots of V^o^'l') is —k, if a = ^ and 



-12a2 + 13a - 2 - ^/72a^ - 168a3 + 123a2 - 28a + 2 
a(2a- l)(6a- 1) 



if I < a 7^ I- But this last expression has values in, say, (—1,0] for | < a 7^ |, implying, by Theorem |l-ll| 
with ^ = 5, that R{^|Ja) < 1 for | < a. 

We focus on the remaining a e (O, ^) \ { '^ g^ } parameter set now. If < a < then ^"(•) has a root 

of the form 

5a - 1 + V-48o3 + 57a2 _ 14a + 1 
2a(2a- l)(6a- 1) 

in the interval (—2 — \/8, —1), meaning that i?(i/'a) < 2 + \/8 here. On the other hand, if < a < | (the 

exceptional value can safely be added again), then V'a(') ^ '^oo* ^-^i^ form 

4a - 1 + \/-8a2 + 8a - 1 
(2a- l)(6a- 1) 

in the interval (— 2 — \/8, — 1 — "\/7), so R{4'a) < 2 + \/8, too, for these a values. (The a=\ case is a 
removable singularity of the above root expression with value —3.) 

Finally, if a = ^^x^, then V'a (-2 - \/8) = 1 - ^ > 0, and, by using (6.1), for any fc > 1 we have 

k 



4 - , then V'a 



-2-^/8)^^^-72^ fc!(A:-l)(fc-2)>0, 



showing that ipa is absolutely monotonic at x — —2 — \/8. But then Theorem 1.12 (applicable because there 
are no poles in (— cx),0]) guarantees absolute monotonicity on the whole interval [— 2 — -^/S, 0], since now 
B{ipa) = +00 trivially, completing the proof of i?3/3,3 = 2 + \/8- 

6.3. Determination of 2 



Formula (2.2) and Lemma 2.8 imply that now we can restrict our attention to rational functions of the 
form 

cz^ + (i - 3a + 3a2) + (1 - 3a) z + 1 
(1 — oz)-' 

with suitable parameters a > and c G M. For simplicity, we have used and will use the letter c instead of 
03. We are going to show that 

V'a,c(-6)>0, V' (-6)>0, C.c(-6)>0 
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and 



imply that a — ^ and c = 2Tt 




Figure 8. The 113/3 2 class. Depicting the roots of 'ipi'^J (•) for fc = 0, 1, 2 as a function of (a, c), the optimal 
point is obtained as the intersection of these three surfaces at (a, c, x) = (g, 2T6'~^)' 



First, by using (2.4 1 evaluated at x = —6, we see that these inequalities are equivalent to 

lOSa^ - 90a - 216c + 13 > 
108a^ - lOSa^ + 42a + 108c - 5 > 
108a'* - 126a^ + 72a^ ~ 12a + 1/2 + (108a - 18)c > 
(a^ - 2a^ + a/2 + c) + fc (-36a'' + 57a^ - 80^ - 36ac - a/2 - 3c) + 

36ac + 2c + 216a^ - 180a* + 26a^ + 216a'^c > 0, 



(6.2) 
(6.3) 
(6.4) 

(6.5) 



where, of course, (6.5) should hold for all fc > 3 integers. We will also use the necessary condition (2.3), which 
now reads as 

-2a^ + a/2 + c>0 (6.6) 



(being just the leading fc-coefhcient of (6.5)) 



Let us express the linear parameter, c from (6.2), (6.3) and (6.6), then combine the resulting inequalities 
to eliminate c and get 

-lOSa^ + 108a2 - 42a + 5 ^ 108a2 - 90a + 13 ^ -2a^ + ia^ - a 



108 - 216 - 2 

Rearranging both inequalities to 0, the two cubic polynomials can be nicely factorized (containing factors 
2a — 1, 6a — 1, 6a + 1 and 36a^ — 60a + 13), so this — together with a > — implies that 



^1 ^5+2^ 

0<a<- or a> . 

6 6 



(6.7) 
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Now let us consider (6.4) and separate two cases according to the coefficient of c. 



If it vanishes, that is, if a = then (6.2) and (6.3) yield c — 2X5, and we are done, because we verify from 



the definition that ?/'a.c(^) — {^^lla is absolutely monotonic on [—6,0] 



Hence the proof is finished as soon as we have shown that (6.7) and a 7^ g l ead to a contradiction. 
First suppose that < a < ^. Then an elementar y ca lculation shows that (6.4 1 (whose left-hand side is 
just (6a — 1) (360"^ — 36a^ + 18a — 1 + 36c) /2) and (6.3) contradict to each off 



As a result, it is sufficient to exclude a > Let us express c (again) from (6.6) and (6.2) to have 

108a2 - 90a + 13 



-2a'^ + 4a^ 



< c < 



216 



(6.8) 



and let us abbreviate the left-hand side of (6.5) by A(fc, a,c) 



Lemma 6.2. For each a > ^+'^"^ and each c satisfying { 6.8), there is an integer k > 3 such that A(fc, a, c) < 



0. 



Proof. Let us pick and fix throughout the proof an arbitrary a > ^+'^"^ and c satisfying (6.8). Since 
dcX{k, a, c) — — (36a + 3)k + 216a^ + 36a + 2, we see that dcX{k, a, c) < holds if and only if 

e32i(a) < fc < e+i(a), (6.9) 



where ^321(0) := |(12a -I- 1) ± ^\/432a2 + 72a + 1. It is also easily seen that — for the allowed a values — (6.9) 
automatically implies fc > 3. So, if k satisfies (6.9), then due to (6.8) we have that 

-2a3 + 4a2 - a\ ^ ^^^^ ^ _ ^a)^ + 35^2 - 8a - l) 



A(fc, a, c) < A yk, a 

Since a > | now, this very last factor is negative if and only if 

36a2 - 8a - 1 



k > 



2a- 1 



(6.10) 



-,f?+i(a) nN 



These mean that (6.9 1 and (6.10) imply A(A:,a,c) < 0. But £'321 (a) < 2a-i~ holds precisely if a e 
(0, 1/6)U (l/2,-|-oo), which is now true by the assumption on a, so to finish the proof, it is enough to 
show that the set 

/36a2-8a-l 
V 2a- 1 

is not empty: an elementary computation yields that for a > ^^g^ we have 

, , , 36a2 - 8a - 1 „ 
f32i(«) 2a-l 

□ 



Remark 6.3. We add that this last integer in the proof, 9, cou7d not be replaced by, say, 10 on the 
whole interval a € (^ ^^^'^ , +ooj . 



Remark 6.4. It is possible to finish the proof Lemma 6.2 by using the dcX{k, a, c) > case, but this 
would result in more difficult computations and would provide (asymptotically) an interval of length 6 
(instead of 9) in the last step. 



7. Determination of R. 



•4/4,2 



We will apply the same principles as in Section 6.3 (or as in the second proof of Lemma 6.1), but, this 
time, we could eliminate the " linear" parameters d and c — recursively via monotonicity — only at the cost of 
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a little bit more computation. 



By appealing to formula (2.2) and Lemma 2.8 again, it is sufRcient to consider rational functions of the 
form 

dz" + cz^ + (i - 4a + 6a2) + (1 - 4a) z + 1 
i^aMz) = 

with suitable parameters a > and c, d G M. (Again, for easier readability we have used and will use c instead 
of as, and d instead of 04.) We show that 

V'a.M(-8)>0, <e,d(-8)>0, VC<i(-8)>0, V'll,(-8)>0 

and 

VfceN,fc>4: V'i^L(-8)>0 

imply a= \, c = ^ and d= 

Let us use (2.4) at a: = —8 (dropping the always positive factors again) and write out the above inequalities 
in detail to get 

Ba + 1) (1 - 4a - 16 ( 6a^ - 4a + - ) + 192c - 2048d ) 



1 - 8(1 - 4a) + 64 ( 6a2 - 4a + ^ ) - 512c + 4096ci > (7.1) 



4a ^1 - 8(1 - 4a) + 64 ^6a^ - 4a + - 512c + 4096(i^ > (7.2) 

lOa^ ^1 - 8(1 - 4a) + 64 (^6a^ - 4a + - 512c + 4096d^ + 
4a(8a + 1) ^1 - 4a - 16 (^6a^ - 4a + + 192c - 2048d^ + 

^ (8a + 1)M 2 ( 6a2 - 4a + M - 48c + 768d] > (7.3) 



10a2(8a + 1) (1 - 4a - 16 ( 6a^ - 4a + - ) + 192c - 2048d] 



2a(8a + ly 1^2 (^6a^ - Aa + - j - 48c + 768dj + -(8a + l)^(6c - 192<i) + 

20a^ ^1 - 8(1 - 4a) + 64 (^6a^ - 4a + - 512c + 4096d^ > (7.4) 

(60" - 6a^ + + 2ac + 2d) + 
(-384a^ + 324a'' - 42a^ - lUa^c - 6ac - 192ad - 12<i) + 
k (4608a^ - 3072a5 + 618a'' + 2304a3c - 36a^ + 144a^c + 4608a2d - 

a^ + 4ac + 576ad + 22d) + 
46080*^ - 2688a^ - 61Ua*c + 300a* - 24576a3d - 4608a2d - 384ad - 12d > 0, (7.5) 



with (7.51 valid for any integer fc > 4. The necessary condition (2.3 1 now reads as 

6a'* - 6a^ + a^ + 2ac + 2d > 0. (7.6) 

The domain of the "non-linear" parameter a > will be divided at the following "natural" points 
(determined by the requirement that denominators in the proof below have constant sign on the corresponding 
intervals) 

„ 3-%/5 1 1 3 3 + %/5 9 + 2^6 

« < < 12 < 8 < 16 < < < ■■■ 
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and also at the "artificial" point . . . < ^ (introduced for technical reasons). The explicit formulae for the 
P42n and g42n expressions (n = 1, 2, . . .) app earing soon are listed in the Appendix (or directly in the proofs). 
The interval a G (^0, '^^^ ^ ■ From (7.4), d is expressed as d < 2048ay-7'6'8a+32 ' "^^ile from (7.3) as d> 
-P422(a,c) -gy joining these inequalities, d is eliminated and we get 



4096a-768 ' 



C < 



-6144a5 + 51200* - 21440^ + 4000^ - 32a + 1 



(7.7) 



From (7.2) we obtain d < ^^^fg^ 
get 



1536a2 - 256a + 24 
which, together with d > ^l^Qa-res ^-^lows us to eliminate d again and 



c > — (-768a3 + 512a2 - 144a + 13) 



(7.1 



From (7.71 and (7.8) we get (8a - l)(8a + l)(16a - 3) < 0, which is impossible for a e [o, ^fi)- 

From (7.4), c can be expressed as c < (2\/5 — 3). But (7.8 1 now says 



1 

384 



The case when a 

1 + 6-\/5) < c, a contradiction 
The interval a G ^ 

^ > P422(a,c) 



16 

, d > 2048a^-76'8a+32 ' (from (7.1 )) d> ^-^^^^^ ■ Moreovcr, an upper estimate d < '^^2048''^^ • 



(containing the unique solution). We have now 3 lower estimates: 



4096a-768 



By making 3 appropriate pairs from these, we eliminate d in each case and get (7.8), 

1 



104a + 11) , 



and 



c < 
c > 



-6144a''+4608a^-1600a''+200a-7 

1536a-128 
-6144a'*+4608a^-1600a''+200a-7 

1536a-128 



if 
if 



16 

i < a < ^ 

12 ^ " ^ 16 ■ 



(7.9) 



(7.10) 



(It is easily seen that a — 1/1 2 in imedi ately leads to a contradiction.) 

On (3 - V5) < a < ^ (7.8) and (7.101 imply |(8a - l)(8a + 1) > 0, which is impossible. 
If < a < ^, then and (| 7.10 [) imply 



c > 
c > 



-61440^+46080-^-16000^+2000-7 



192 



15360-128 

-768a3 + 512a2 



144a + 13) 



if 
if 



12 < O. ^ 8 ' 

i < a < ^ 

8 ^ " 16- 



(7.11) 



On I < a < (7.9) and (7.11 ) cannot be true simultaneously. On the other hand, if < a < |, then (7.9) 
and \7.n\ yield a = 1/8 and c = 1/128, and with d > ^^^^ and d < ^'^a^"^ we get d = 1/4096. 



The interval a G 



_3_ 9+2^6 \ 
16' 24 J' 



We will use here the two lower estimates d > ^^^q''^^''' - d > 



-Ga'^ + 6a-* 



2ac) (this second one derived from (7.6|) and the upper estimate d < 



P424(a,c) 

4096 

P423(g,c) 



2048 



. Again, 



by forming two appropr iate pairs, d is eliminated and we obtain two inequalities with a and c. After a 
rearrangement, we have (7.9 1 and c > (—7680'^ + 768a-^ — 160a + 7), leading to a contradiction. 



So far we have examined the region < a < ^^24^ ^^'^ proved that the system has a solution if and only 
if a= 1/8. 

In the rest of this secti o n, we will show that there are no solutions if a > ^'^24^ ■ '^^ '"^-'^ shown (in a few 
lines, by using only (7.1), (7.2) and (7.6), but skipping the details here) that if a > ^^^/^ , then we have 




768a2 
-a2 



160a + 7) < c < 



12 



6a3 - a^ - 2ac) <d< ^ {768a^ 



24 

-1536a3 + 1728a2 - 424a + 25) , 
- 512a2 - 512ac + 104a + 192c - 7) 



(7.12) 



or 




-1536a3 + 1728a2 - 424a + 25) < c < (l92a2 - 104a + ll) , 
ike (-384a2 + 224a + 512c - 25) < d < ^ {768a^ - 5120^ - 512ac + 104a + 192c - 7) 



(7.13) 



Let us introduce the abbreviation A(fc, a,c, d) to denote the left-hand side of (7.5). 
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Lemma 7.1. For each ^(9 + 2^/6^ < a < ff, and any c and d satisfying {7.12], we have that 
A(77, a,c, d) < 0. 



Proof. Since the polynomial ddX{77,a,c,d) — — 48^425(0) has a unique root £1421 ~ 1.17854 in the given 
a-interval, we separate two cases. 



The interval a G (9 + 2v^) , £1421) . Here 5dA(77, a, c, d) > 0, so by ( 7.12 ) 



A(77, a,c,d)<X[ 77, a, c, {768a^ - 512a^ - 512ac + 104a + 192c - 7] 



57 
128* 



3a + l)P426(a,c), 



from which we see that it is enough to show that P 426(a, c) < 0. But for the current a and c values we have 
dcP426{a,c) = 192 (I28a2 - 800a + 925) > 0, so by ^LU) 



P426(a, c) < P426 ( a, — (-1536a3 + 1728a2 - 424a + 25^ 



;(8a + 1) (24576a'* - 139264a^ + 2122880^ - 126000a + 17575) < 0. 



The interval a e [£'421, f§]. Here ddX{77,a,c,d) < 0, so by (7.12) 

A(77, a, c,d) < X (^77, a, c, ^ (-60"* + 6a^ - - 2ac) j = 12a(8a + l)P427(a, c), 

so, again, it suffices to show P427(a, c) < 0. But now we observe that dcP427{a,, c) — 192a ^ — 18 24a + 2850 is 
positive, if £1421 < a < | (38 - \/494), and non-positive, if | (38 - ^/494) < a < f§. By (7.121 again, in the 
first case, we have 



P427(a, c) < P427 ( a, — (-1536a^ + 1728a2 - 424a + 25) | = 



1 
256 

and in the second case 



(8a + 1) (61440^ - 629760^ + 1744960^ - 172672a + 35625) < 0, 



^'427(0, c) < P427 ( a, ^ (-768a3 + 768a2 - 160a + 7] 



32^ 



ia + 1) (256a2 - 648a + 175) < 0. 



□ 



Remark 7.2. Mathematica tells us that among the integers 4 < fc < 1000, exactly members of the 
interval 77 < fc < 98 share the property that X{k,a,c,d) < for any allowed a, c and d in the previous 
lemma. 



Lemma 7.3. For each ^ (9 + 2\/6) !i ^ f§i and any c and d satisfying \7.13V, we have that 
A(54, a,c, d) < 0. 



Proof The proof is completely analogous to the previous one. Now (9dA(54, a, c, d) — — 24P428(a) has two 
roots ^p422 ~ 0.808208 and P423 ~ 1.97947 in the given a-interval, so we again separate two c ases. 
The interval a € [54 (9 + 2V&) , Q422) U (£(423, i§] • Now 9dA(54, a, c, d) > 0, so by (|7.13) 



A(54, a, c, d) < A ^54, a, c, (768a^ - 5120^ - 512ac + 104a -I- 192c - 7) 
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159 
256 



(8a+ l)P429(a,c). 



We sh ow th at P429 < 0. Indeed, 5cP429(a, c) = 64 (l92a2 - 832a + 663) is positive if a e (9 + 2^/6) , £1422) 
so by (7.13) here 

P429(a,c) < P429 (^a, ^ (I92a2 - 104a + U] 
hsa + 1) (12288a'' - 58368a3 + 761600^ - 35152a + 4199) < 0. 



On the other hand, dcPi29{a, c) < 0, if a e (^423, f§] , hence by (7.13) we have now 
-P429(a, c) < P429 (a, ^ (-1536a^ + 1728a2 - 424a + 25) 

^(8a + 1) (12288a'' - 46080a^ + 53888a2 - 29328a + 4199) < 0. 



The interval a G [9^22, 9^23]- Here 9dA(54, a, c, d) < 0, so by (7.13) 

A(54,a,c,d) < A (^54, a, c, ^ ( 



384a^ + 224a + 512c - 25 = 



512 



(8a + l)P42io(a,c). 



By examining the sign of 9cP42io(a,c) = -512 (256a3 - 5088a2 + 16536a - 11713), we prove finally that 
P4210 < as well. We see that 9cP42io(a,c) > 0, if ^422 < a < £1424 ~ 1.00126, but 9cP42io(a,c) < 0, if ^424 < 
o- < f?423- By taking into account ( |7.13 ), in the first case we have 

P42io(a, c) < P4210 [a, ^ (I92a2 - 104a + ll)^ = 

53(8a + 1) (I2288a4 - 58368a^ + 76160a2 - 35152a + 4199) < 0, 
while in the second case 

P42io(a, c) < P4210 {a, ^ (-1536a3 + 17280^ - 424a + 25) 

8a(8a + 1) (6144a'* - 45312a^ + 102368a2 - 86284a + 17225) < 0. 

□ 



Remark 7.4. Accoiding to Mathematica, fc = 54 is the only integer in the interval [4, 1000] H N such 
that A(fc, a, c,d)<0 holds for any admissible a, c and d triples in the previous lemma. 



Lemma 7.5. For each || < a, and any c and d satisfying ( 7.12), there exists an integer k > A such that 



A(fc, a, c, d) < 0. 



Proof. Let us fix any admissible a, c and d throughout the proof. We first compute ddX{k,a,c,d) = 
2P42ii(fc7a)i where 

P42ii{k, a) = k^ - (96a + 6)k'^ + (23040^ + 288a + ll) fc - 12288a3 - 2304a2 - 192a - 6, (7.14) 

and show that this cubic polynomial (in fc) has three distinct real roots and also give some bounds on the 
roots in terms of a. 

It is seen that the function fc 1— >■ 9fcP42ii(fc, a) has two distinct real roots at fc = ^425(1) ■— 
I (96a + 6 ± \/3v'768a2 + 96a + l) , being the abscissae of the two strict local extrema of fc i~> P42ii(fc,a). It is 
also easily established that P42ii(P425(a)j «) > £^nd P42ii(f?425(«)j o-) < 0- By observing that P42ii(4, a) < 
and 4 < ^425 (a), further, by taking into account that the leading coefficient of fc 1— > P42ii(fc, a) is positive, we 
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get that each of the three disjoint intervals 

(4, ^?425(a)) U (fi-rasCa)' (?425(a)) U (^?i)5(a)' +00) (7-15) 

contains precisely one root of fc H' -P42ii(fc,a) = ^9dA(fc, a, c, c?). 

By denoting these roots by £'426(11)1 < 042e{O')2 < £'426(0)31 we also see from the properties of the cubic 
polynomial that i f k € (g426(Q)2; £"426(0)3) 1 then dd^{k,a,c,d) < 0. (Note that for any such k, k > A holds.) 
Also referring to (7.12), this means that if fc G (£'426(0)2, £'426(0)3), then 

A(fc, a, c,d) < X ^fc, a, c, - (—6a* + 60"^ — — 2ac)^ = 6a(8a + 1)^4212(^1 o, c), 

where 

P42i2(/c, o, c) = fc^ (4a3 - + a + c) + fc (-192a* + 212a^ - 270^ - 48ac - 2a - 3c) + 

1536a^ - 1344a* + 104a^ + 384a2c + ISa^ + 48ac + a + 2c. 

Clearly, in order to finish the proof, it is sufficient to show that there is an integer k g (5426(0)27 £'426(0)3) 
such that /42i2(fc, o,c) < 0. 

Let the two roots of /c dcP42i2{k, o, c) = fc^ - (48a + 3)k + 384a^ + 48a + 2 be denoted by ^427(0) := 
i (48a + 3 ± V768a2 + 96a + l) , then dcPi2i2{k, a, c) < if A: S (£'^27(0), £1427(0)) • We wih show next that 

(5426(0)2, £'426(0)3) n (5427(0), £'ii7(o)) = (5426(0)2, 5i27(a))- (7-16) 

In fact, we have seen in ( |7.15[ ) that 5425(0) < 5426(0)2, said an elementary calculation shows that 54^27(0) < 
5425(0)- An analogous argument tells us that 5^27(0) < 5425(0) < 5426(0)3, verifyin g ( 7.16). 

Thus, by picking any k e (5426(0)2, 5427(0)), we can conclude (with the help of (7.12) also) that 



Pi2i2{k, o, c) < P4212 (^'^'^92 (^'^'^^"'^ 



768a2 - 160a + 7 



192^ 



Ba + l)(fc - 1) ((7 - 24a)A: + IQSa^ - 96a - 14) 



from which we see that the proof is finished if we find any k G (5425(0)2, 5427(0)) H N such that (7 — 24a)fc 
768a2 - 96a - 14 < 0, or, in other words, k > 

(Notice that this time both 5425(0)2 < 5425(0) and 



24a-7 

We aim to show now that 5425(0)2 < '^^^"'^ia-r'^'^ ■ 
''^^°'2Aa-7^^^ < 5^^25(0) are true, so we have to take an extra step.) But by referring to (7.15) we have 
5425(0)1 < 5425(0), and an elementary computation yields that 5425(0) < ^^^^ 
hard to see that 



24a-7 



-, further, it is not 



4211 



/768a2 - 96a- 14 
V 24a- 7 



,0 = 



48a(8a + l)(16a - 3) (9216a3 - 1512a + 49) 
(24a - 7)3 



< 0. 
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However, we already know that P42ii(-, a) < if and only if fc G ( — oo, g426{a)i) U (£'426(0)2, £'426(0)3)1 so 
these force ^^426(a)2 < "'^Va"-'?"'^ ■ 

Summarizing, it is enough to prove as a last step that =/= ( -^^^^^j^^^f— ^ , f?427(o)) being true, since 



Remark 7.6. Similarly to Remark 6.3 the bound 9 could not be replaced by, say, 10 on a £ (22/10, +00). 



Remark 7.7. The proof we have just presented for Lemma 7.5 (as well as the one we will present for 
le next lemma) would break down on the larger in 
finer subdivision a e (9 + 2\/6) , f|] U (f|, +00) 



the next lemma) would break down on the larger interval a > ^^2i^ ■ This explains why we have chosen the 



Lemma 7.8. For each < a, and any c and d satisfying {7.13), there is an integer fc > 4 such that 
A(fc, a, c, d) < 0. 



Proof We have seen in the preceding proof that ddX{k, a, c, d) > if fc € [£?426(a)i, £'426 (a)2] C (4, +00), 
so — by (7.13) — for these fc values we have A(fc, a, c, d) < 



A fc, a, c, 



2048 



(768a3 - 512a^ - 512ac + 104a + 192c - 7) = 



Ba + l)(fc - 1) 
1024 



P42i3{k,a,c) 



where 



P42i3(fc, a, c) = fc2 (768a3 - 7680^ + 160a + 192c - 7) + 
fc (-49152a^ + 39168a3 - 5376a2 - 12288ac - 128a - 960c + 35) + 
589824a5 - 294912a'' + 291840^ + 147456a2c + 3840a2 + 24576ac - 384a + 1152c - 42. 
To show that P42i3(fc, a, c) < for a suitable fc G [£1426(0)11 £'426(0)2] H N, we first compute 
9c^42i3(fc, a, c) = 192 (fc2 - (64a + 5)fc + 768a2 + 128a + 6) , 

and denote the roots of fc 1— !■ 9c^42i3(^, o, c) by g,± g(a) := i (64a + 5± Vl02 4a^ + 128a + l). Note that 
dcP42i3{k,a,c) < if £'4^28(0) < fc < £'428(0). Now by using formulae (7.151 and (7.14), one easily shows that 
£'425(0) < £'428 (a) < £'425 («) and 



-P4211 (£1428(0), o) = 32a(8a + 1) ( Vl024a2 + 128a + 1 - 16a - l) > 0, 



implying £'426(0)1 < ^428(0) < £'426(0)2- Hence if fc e (g^^g^ia), £'426(0)2] n N, then, also by ( 7.13 ) 

1 



-P42i3(fc, o,c) < P4213 k,a, 



512 



-1536a^ + 1728a^ - 424a + 25) ) = i(8a + l)P42i4(fc, o 



with 



P42i4(fc, a) = fc^ (I92a^ - 144a + 19) + fc (-12288a^ + 21120^ + 1296a - 95) + 



147456a'' + 110592a^ - 27264a^ - 2016a + 114. 

Now the leading coefficient, 192a2 — 144a + 19 of fc !-> P42i4(fc, o) is positive, and its discriminant, A42(a) := 

37748736a^ - 51904512a5 + 46043136a'* - 14751744a^ + 2101632a2 - 27360a + 361 

is also positive. (This last statement can easily be shown by noticing that A42 (o) > for all a £ M, and 
A42(22/10) > 0, so A42(a) > for a > 22/10. Repeating this recursively and similarly for the lower order 
derivatives we get that A42(a) > for all a > 22/10.) This means that we can denote by ^^29(0) the two real 
roots of the quadratic polynomial fc 1— >■ P42i4(fc, o). Clearly, £'429(0) < fc < ^429(0) implies P42i4(fc,a) < 0. 
In light of the above, in order to finish the proof it is enough to show that 



<^ ^ (£'429(0), £'429(0)) n (^£'428(0), ^426(0) 



nN. 
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More specifically, we can prove that the following sufficient condition 

Q42ai^) < 6'429(«) < 16a + 19 < 16a + 24 < £1426(0)2 < £'429(0) (7-17) 
holds: the first inequality is true because of 

^'4214 (£'l28(«)' a) = 128a ((24a - 7)v/l024a2 + 128a + 1 + 136a + > 

and 

dkP4214 (f?428(o)iO) = -6144a2 + 1792a - (l92a2 - 144a + 19) v/l024a2 + 128a + 1 < 0; 

the second inequality is a consequence of 

P42i4(16a + 19, a) = -16 {576a^ - 12640^ + 1209a - 323) < 0; 

the fourth inequality follows from P42ii(16a + 24, a) = 40960^ + TeSa^ - 25360a + 10626 > and 16a + 24 < 
£'425(0); while the fifth inequality is derived from the fact that 

+ _ 12288a3 - 2112a2 - 1296a + 95 + \/^42(o) 
^425^"^ < ^^429(«) - 2(192a2-144a+19) 

and 

r> / + r \ \ 96o 

P4211 £'490(0)1 o = 5- X 

^^''^ ' ' (192a2-144a+19)' 

(^301989888a*^ - 622854144a^ + 474218496a6 - 1089110016a5 + 816156672a'* - 250172928a^+ 

31469760a2 - 960184a + 9025+ 

(98304a^ - 1351680^^ - 132096a^ + 106880a2 - 19464a + 475) ^/A^a)J < 0. 

□ 

Remark 7.9. In the proof of this lemma — in order to establish A(fc, a, c, d) < — one would have 4 initial 
possibilities to choose from according to the signs of the appropriate da and dc partial derivatives. It turns 
out that if one chooses certain 2 out of these 4 possibilities, then the proof could not be completed. On the 
other hand, the remaining two feasible directions (out of which we have presented one) require approximately 
the same amount and type of computations. 



Remark 7.10. The motivation for {7.17) came from the fact that 



hm ^429(0)^ 16 < lini ^^2^ =root2(l, -96, 2304, -12288)^ 26.4432. 

a— !.+oo a a— S.+00 a 

In fact, for a > 22/10, we have £'426(0)2 — £'429(0) > 6, but . . . > 7 does not hold. 



8. Determination of RA/i,p for 3 < p < 5 

For the sake of brevity and due to the fact that the techniques here are very similar to the corresponding 
earlier ones, we provide only brief (indication of the) proofs of the remaining results. 

8.1. Determination of i?4/4 5 

We prove here that -R4/4_5 ~ 3.7432. The exposition is completely analogous to the one in Section 
There are exactly 4 rational functions in this 114/4 5 class: 

a^z^ + a^z^ + a2Z'^ + aiz + 1 
(1 — az) 



6.1 



^(^) = n „„,4 
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where 

a4 =root,„,(2160000,4032000, -93600,-120,1), = root™, (18000, 24000, -1800, 0, 1), 
02 = root™,(1200, -7200, 120,200,3), ai = root„,^ (15, 60, -30, -20, 7), 
a = root,„„(120, -240, 120, -20, 1), 

with 

(m4,m3,m2,mi,mo) G {(1, 1, 4, 1, 4), (2, 2, 2, 3, 2), (3, 3, 3, 4, 1), (4, 4, 1, 2, 3)}. 

First we remark that in aU 4 cases the only pole of ip is positive real. 

If 7714 = 1, then the corresponding tp satisfies ip^^'>{— 1/2) < < ^^^\0) — 1, hence t/;^^^ has a root in 
(-1/2,0), so by Theorem [m] with ^ = 5 it has R{ip) < 1/2. 

If m4 = 2, then the corresponding -0 satisfies -tjA^^O) < < ^'^^•'(^1): hence by Theorem 1 1 . 1 1 1 with £ = 8 
it has R{ip) < 1. 

If TO4 = 4, then the corresponding ip satisfies < < ?A'^-'(0), hence by Theorem 1 1 ■ 1 1 1 with £ = 7 

it has R{-ip) < 1/2. 

If TO4 = 3, then the corresponding ip satisfies ^'(x*) — with 

a;* := rooti (1, 12, 48, -32, -864, -2016, 2784, 13248, -9072, -35136, 44928, -20736, 3456) w 

-3.743299247417768882803493, 



hence by Theorem 1.11 with ^ = 1, it has R{^) < |a;*|. By using partial fraction decomposition, we can show 
that R{ip) = \x*\, therefore i?4/4,5 = 

8.2. Determination of R^/i^i 



We prove that i?4/4.4 ~ 5.1672. We refer to Section 6.2 for an analogous treatment and for further details 



on the theorems used. Elements of this class can be written as 

^a{z) = (1 - az)-^ [a^ - 4a3 + Sa^ - ^ + ^) + (^-4a^ + Qa" - 2a + i^) 

6a2 - 4a + + (1 - 4a)z + 1 

with a real parameter a. Here also we can assume that a > and ipa does not have removable singularities 
(indeed, a = implies i? = 1; a < without removable singularities implies i? = 0; while if the numerator 
and denominator have a common root, then ipa is one of the 3 rational functions in Section |6.1[ hence 
R ^ -^3/3,4 < 5). By applying the same techniques as earlier in the one-parameter families, we can prove 
that 

^4/4 4 = root3(l, -24, 240, -1168, 1848, 7008, -30528, 7488, 71568, 36864) w 

5.167265421277419938673374, 
with the unique optimal rational function corresponding to the parameter value 

a = a* := rooti (147456, -546624, 799488, -601344, 258432, -66576, 10400, -960, 48, -1) w 

0.09713312764144710280835106. 

8.3. Determination of i?4/4^3 

This class ca n be described by 2 real parameters. We can prove along the lines of the uniqueness-type 



proof in Section 6.3 that i?4/4,3 = 3 -|- VlS « 6.8729. 
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9. Appendix 

In this section we list some auxiliary algebraic numbers, polynomials or functions that have appeared in 
certain more involved proofs. 

9.1. Algebraic expressions in the proof of R4/4J 

(?47i =rooti(252105,936390,1441629,1175608,534576, 128352, 12704) w -0.843194, 
(3472 = root2(252105, 936390, 1441629, 1175608, 534576, 128352, 12704) w -0.471357, 
(3473 = root2(l, -30, 390, -2760, 11160, -25200, 25200) w 7.64527, 
(?474 = rooti(l, -15, 90, -210) w 5.64849, 

P47i(a) = 4593387934777490821322994ai° - 234741768165039984427600980^+ 
54351276637181597088697031a* - 75708182541462463946985360a^+ 
717915443131527432118064640^ - 513392369081631351916955400^+ 
323834997962359967667069780^^ - 22447171953065668650772896o3- 
31555722050640962925153690o2 - 13278623731050880370830074o- 
1667865061502294482628289, 
£1475 = rooti (2914539265575, 14876524399779,33497711997246,43683499824678, 
36368607954483, 20052618149655, 7324770907832, 1709817444048, 
231520801344, 13859993824) w -0.850052, 
(?476 = root2(2, 0, -570, -6360, -28755, -58950, -44775) w 21.5907, 
(3477(a) = rooti(343o3 + 4410^ + 189o + 27, -61740^ - 76440^ - 3150a - 432, 

473340^ + 55860a2 + 21882a + 2844, -1975680^ - 2199120^ - 80976a- 
9856, 4784850^ + 5005350^ + 173460o + 19950, -6482700^ - 6438600^- 
215460a - 24480, 3858750^ + 3682350^ + 120960o + 13440), 
£1478 = rooti (1525366344600, 8651006268660, 22401405617094, 35008454688795, 
36765395913141, 27330042736584, 14744440158000, 5816571277965, 
1665232178673, 337417309264, 45932541984, 3771903744, 141310592) « 
-0.469514, 

£1479 = root3(1715, 2793, 1428, 232) w -0.358565, 

£14710 = rooti (3675, 3507, 1152,128) w -0.274796, 

(34711 (o) = rooti (7o + 3, -84a - 32, 420o + 120, -840a, -840), 

£-4712(0) = root2(7a + 3, -84a - 32, 420o + 120, -840o, -840), 

a*7 = rooti (3690134492416632557532940385381064525, 
61029736485482612404619408048530290750,486403976093629932677089384897095046935 
2487196669730164018264503869937911017740,9165142322735094610456262590295894934837 
25915196147894075773838728122892827040386,58458916104726476816514027882738196455735, 
107979092136285121513038059062379021459880,166321122780685405075074427832862986958975, 
216449429305753337078085597912060934057890,240240731812681651181007922270907964916605, 
228911785323865533798715657213902093130380,188040840503767270695413477386184055716415, 
133446371333215106255009916213506161760670,81810565480540039001637010346621868861045, 
43219005242542299485983226045186624870400,19562352897200734179995371112351361862080, 
7506336744004524872039975974458313909440,2394351571156175328455192580869480281280, 
610166214096532168321831058806216604160,112156289536490108521784144024335180800, 
8974819895289125868510161871145369600,-2895655868721204798705192374322355200, 
-1556648822533171464739057080454594560,-420567479322606922282373395209093120, 
-79380339166868822139598939104215040,-11057440115321274264887170064056320, 
-1128821021046790224548013859143680,-80578528495036917551879568752640, 
-3612945212567877804734004854784, -76850622516036469593789693952) « -0.439849, 

B (0^7) = root2(4191472, 370695456, 15701398968, 423572490288, 8166117227943, 
119697694352106, 1385540391042992, 12982888147808790, 100093600309232610, 641253042735937920, 
3429070908298155495, 15292307228951973150, 56488604555080263600, 170258230386661619700, 
405691319555093173950,698620766882164002000,475967180133964116000, 
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-2768737485607368840000,-19171364099094461844000,-83177899383693915360000, 
-273285810570616506720000, -658873655023431060000000, - 1076077292526138186000000, 
-1039429806316804224000000,-402385174364630880000000,1632495959008528320000000, 
15171036284831515200000000,63475278812225280000000000,141035850392839718400000000, 
165407760061818624000000000, 81912466393111872000000000) w 2.743911. 

9.2. Algebraic expressions in the proof of i?4/4.2 

P42i{a, c) = 1536a5 - 1280a-* - 5120^0 + 4360^ + 5760^0 - 54a'^ - 72ac + 2a + c, 

P422(a, c) = 2304a'* - 1728a3 - 10240^0 + 480a2 + 768ac - 48a - 48c + 1, 

^'423 (a, c) ^ 768a3 - 5120^ - 512ac + 104a + 192c - 7, 

-P424(a, c) = -384a2 + 224a + 512c - 25, 

^425 (a) = 512a3 - 7296a2 + 22800a - 17575, 

(3421 = rooti(512, -7296, 22800, -17575) w 1.17854, 

^426 (a, c) = 983040^ - 6799360* + 1266432a3 + 24576a2c - 827264a2 - 153600ac + 156400a+ 
177600c - 6475, 

f427(a, c) = 768a5 - 8064a* + 200720^ + 192a2c - 158530^ - 1824ac + 2888a + 2850c, 

P428(a) = 1024a3 - 101760^ + 22048a - 11713, 

£1422 = rooti(1024, -10176, 22048, -11713) « 0.808208, 

£1423 = root2(1024, -10176, 22048, -11713) « 1.97947, 

^429 (a, c) = 49152a5 - 245760a* + 3653120^ + 12288a2c - 210496a2 - 53248ac + 38272a+ 
42432c - 1547, 

f42io(a, c) = 5406720a5 - 28110848a'' - 131072a3c + 44025856a3 + 26050560^0 - 269664000^- 

8466432ac + 5517512a + 5997056c - 292825, 
g424 = rooti(256, -5088, 16536, -11713) w 1.00126. 
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